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We consider noisy, non-Iocal unitary operations or interactions, i.e. tlie corresponding evolutions are de- 
scribed by completely positive maps or master equations of Lindblad form. We show that by random local 
operations the completely positive maps can be depolarized to a standard form with a reduced number of pa- 
rameters describing the noise process in such a way that the noiseless (unitary) part of the evolution is not 
altered. A further reduction of the parameters, in many cases even to a single one (i.e. global white noise), is 
possible by tailoring the decoherence process and increasing the amount of noise. We generalize these results 
to the dynamical case where the ideal unitary operation is given by some interaction Hamiltonian. The result- 
ing standard forms may be used to compute lower bounds on channel capacities, to simplify quantum process 
tomography or to derive error thresholds for entanglement purification and quantum computation. 

PACS numbers: 03.67.-a, 03.67.Pp, 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 

Quantum systems evolve via unitary operations U{t), as 
they are governed by the Schrodinger equation. This also 
holds for composite systems, e.g. a small quantum system 
S which is surrounded by some environment E, where the 
evolution of the total system is described by a unitary opera- 
tion Ussit). The dynamics of the system S alone — which 
can be obtained by tracing out the (uncontrollable) degrees of 
freedom of the environment — is in general no longer unitary. 
In fact, the system interacts with degrees of freedom of the 
environment, leading to entanglement between system S and 
environment E reflected in UsE{t) ^ Us{t) ^ UE{t). The 
system-environment interaction leads to decoherence and the 
dynamics of the system can be described either by a (time 
dependent) completely positive map (CPM) £{t) or — under 
certain assumptions on the nature of interaction — by a mas- 
ter equation of Lindblad form |1]. From the perspective of 
quantum information processing, such an interaction with en- 
vironmental degrees of freedom is undesirable and leads to 
errors and noise in the system. As discussed below, an arbi- 
trary noise process acting on a d-dimensional system S is, at 
fixed time t^, determined by 0{d'^) real parameters. Even for 
small system sizes, e.g. when S consists of three qubits (i.e. 
d = 8), this leads to a huge number of independent parameters 
(e.g. around 4000 for the three-qubit system), which makes 
an analytical treatment of the influence of such general noise 
processes on the properties of the system rather difficult. This 
is particularly hindering when considering either large sys- 
tems or sequences of several noisy evolutions (or gates), as is 
e.g. required in the analysis of quantum circuits or processes 
such as entanglement purification. 

When considering the influence of noise in quantum in- 
formation processing, one hence often restricts the analysis 
to certain (ad hoc) noise models, such as Pauli channels or 
depolarizing (white) noise models. This is usually the case 
in the analysis of entanglement purification protocols in the 
presence of noisy operations as well as in the theory of fault- 
tolerant quantum computation. On the other hand, having a 



specific physical set-up in mind, one can sometimes justify 
these (or other) noise models by a microscopic description of 
the underlying system-environment interactions, where only 
the dominant part of noisy interactions is considered. How- 
ever, when considering (abstract) quantum processes that deal 
with the manipulation of quantum information, one does not 
want to restrict oneself to a specific physical set-up, but rather 
would like to keep the analysis at an abstract level and as gen- 
eral as possible. To this aim, it would be very useful to justify 
the usage of simple noise models in a general context, or to 
provide a method to bring any noise process to a simple stan- 
dard form described by a few parameters. 

In this article, we provide a method which allows one to 
achieve this aim. We show that indeed any noise process can 
be brought to a simple standard form by means of depolariza- 
tion. That is, by applying appropriate (local) unitary control 
operations on a system before and after the noisy evolution in 
a correlated way, one can depolarize the noise process. This 
depolarization process of the corresponding CPM £ or the Li- 
ouvillian C can be viewed as an analogue of the depolariza- 
tion of mixed states. For bipartite states, for instance, it was 
shown that one can bring any mixed state p of two d-level sys- 
tems to a standard form specified by a single parameter using 
an appropriate (random) sequence of local unitary operations. 
Depolarization takes place in such a way that the fidelity of 
the state, i.e. the overlap with a maximally entangled state 
1$) — l/\/rfX]fc=i 1^)1''')' remains invariant. The resulting 
states (isotropic states) are equivalent to Werner states and 
are given by p{x) = a;|<I>)($| + (1 — x)-^l. Werner states 
played an important role in the investigation of the relations 
between entanglement and local hidden variable theories 01, 
as well as in the development of entanglement purification 
protocols, schemes which are becoming increasingly impor- 
tant since it was realized that entanglement can serve as a valu- 
able resource not only in quantum communication but also in 
quantum information processing. The development of these 
important issues was triggered by the simplified description 
of Werner states (still covering essential entanglement prop- 
erties), and allowed at the same time to obtain necessary or 
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sufficient conditions for separability or distillability for arbi- 
trary bipartite states. 

We are confident that also the depolarization of noisy evo- 
lutions will prove to be a fruitful tool in the analysis of noise 
processes. In direct analogy to the depolarization of states, 
the depolarization of the noise maps takes place in such a way 
that the fidelity of the ideal (unitary or Hamiltonian) part of 
the evolution is not altered. In fact, we make use of the iso- 
morphism between completely positive maps and mixed states 
fs*, "Tl and connect the problem of depolarization of maps to 
the depolarization of the corresponding states, while respect- 
ing certain locality restrictions. For decoherence processes 
(e.g. storage errors of a system due to its interaction with the 
local environment, or errors resulting due to sending a sys- 
tem through a noisy quantum channel), where the ideal oper- 
ation is the identity or, equivalently, the Hamiltonian part in 
the corresponding master equation is zero, we find that one 
can depolarize the corresponding map or master equation to 
a standard form which is described by correlated and uncor- 
rected white noise processes. In the case of two d-level sys- 
tems, for instance, the corresponding depolarized map is de- 
scribed by three real parameters (the weights of ideal opera- 
tion, single particle (uncorrected) white noise processes and 
two-particle (correlated) white noise) as compared to 0{cfi) 
parameters of an arbitrary map. We also consider noisy inter- 
actions (i.e. the ideal evolution is given by some non-local 
unitary operation U or some non-trivial system interaction 
Hamiltonian H), where we concentrate on two-system inter- 
actions. We find that for certain unitary operations, in par- 
ticular for SWAP gates, CNOT gates, as well as phase gates 
with arbitrary phase a, a depolarization is still possible. The 
required number of parameters to describe the (depolarized) 
noise process depends on the unitary operation (interaction) 
that has to be kept invariant, and is given by 17 in the case of 
arbitrary phase gate, 8 in the case of the CNOT gate and 3 for 
the SWAP gate. 

Knowledge of the exact form of the noise process (which 
can e.g. be acquired by means of gate tomography) or addi- 
tional control of interactions (e.g. the ability to switch a noisy 
interaction on and off at will) allows to further tailor the noise 
process. In this case, the fidelity of the ideal operation is de- 
creased by a certain (small) amount, while the description of 
the noise process is simplified and the number of relevant pa- 
rameters is further reduced. In many relevant cases (e.g. noisy 
SWAP or CNOT gate, switchable noisy phase gate), one finds 
that one can indeed simplify the noise process in such a way 
that the corresponding CPM is described by a single parameter 
and the noise process corresponds to correlated white noise. 
The total amount of noise is — in the worst case — increased 
by about an order of magnitude, as weight of the ideal oper- 
ation is transferred to the noise part in an appropriate way to 
achieve this further simplification. 

While in the case of maps a depolarization with a signif- 
icant reduction of the associated parameters is only possible 
for certain unitary operations, one finds that, in the case of 
master equations, sequential application of fast intermediate 
local unitary control operations allows one to depolarize any 
master equation (of two systems) to a standard form described 



by at most 17 parameters. In return this depolarization pro- 
tocol generally increases the noise level of the decoherence 
process. Under certain circumstances, one may even achieve 
a standard form described by a single parameter for arbitrary 
two-qubit interactions by accepting a further increase of the 
noise level. 

Such standard forms for noisy evolutions may have wide 
spread applications in the analysis of quantum information 
processes under realistic conditions. For instance, our ap- 
proach allows one to obtain lower bounds on the capacity 
of arbitrary multipartite quantum channels by considering the 
corresponding depolarized channels. The depolarized noise 
process also gives rise to a simplified process tomography. 
The tomography has to reveal fewer parameters (the param- 
eters characterizing the standard form) than those necessary 
to describe the original decoherence process or the noisy gate. 
This can lead to a significant reduction of the experimental 
effort to sufficiently characterize the influence of noise in a 
given set-up. Also processes involving sequences of noisy 
gates, e.g. entanglement purification or some quantum cir- 
cuits, can be analyzed by considering the standard forms for 
the corresponding gates. The resulting threshold values do 
no longer refer to specific error models but are valid in gen- 
eral, as any noise process can be brought to the corresponding 
standard form. When applying this method to derive gener- 
ally valid error thresholds, e.g. in the context of fault tol- 
erant quantum computation, some care is required. An im- 
plicit assumption in order to allow the application of such a 
local depolarization procedure is that the corresponding (lo- 
cal) control operations can in fact be (noiselessly) applied to 
the system. When dealing for instance with decoherence pro- 
cesses due to channel noise or local interaction of the system 
with some environment, such an assumption is perfectly rea- 
sonable. Also for two-system interaction gates (such as the 
CNOT), one may assume that local, single system gates are 
noiseless (or introduce a negligible amount of noise as com- 
pared to the two-system gate). However, when dealing also 
with noisy single system operations (as is e.g. required in the 
analysis of fault tolerant quantum computation), it is no longer 
straightforward to apply our results. One might argue that for 
sequences of gates the required (random) operations for de- 
polarization can be incorporated in previous/subsequent noisy 
gates, although it is not entirely clear whether this argument 
justifies the assumption that any gate within a quantum circuit 
is already of standard form. However, whenever local, single 
system control operations can be assumed to be noiseless, our 
results are applicable and one can indeed bring an arbitrary 
(non-local) noise process to a simplified standard form. 

This paper is organized as follows: In Sec.|ll|we review 
basic properties of the Jamiolkowski isomorphism between 
completely positive maps and states, which will be the main 
tool for the derivation of standard forms for CPM in the fol- 
lowing sections. We will then apply this isomorphism in 
Sec. |in] in order to provide standard forms for an arbitrary 
decoherence process in the case where the corresponding con- 
trol operation to achieve this standard form does not have to 
obey any locality requirements. In Sec. lIVI we derive standard 
forms for maps describing arbitrary decoherence processes 
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and some noisy unitary operations. These standard forms are 
achieved by control operations that are local with respect to 
(w.rt) some given partitioning. In Sec.|V]we suggest a proto- 
col to bring an arbitrary noisy evolution described by a master 
equation into some standard form, for which the accompany- 
ing noise process is described by a reduced number of param- 
eters. Finally we summarize our results in Sec. IVll Some 
technicalities can be found in the appendices. 



II. THE JAMIOLKOWSKI CORRESPONDENCE 
BETWEEN COMPLETELY POSITIVE MAPS AND STATES 

In this section we review some properties of the 
Jamiolkowski isomorphism jsj between completely positive 
maps (CPM) and states. In Sec. Ill Al we state and discuss 
this isomorphism first on an abstract level as a correspon- 
dence between matrices and the endomorphisms of the cor- 
responding matrix algebras. We will then restrict this gen- 
eral isomorphism to the physical setting of quantum states 
and quantum operations in Sec. IIIBI where the isomorphism 
has a clear interpretation in terms of a teleportation proto- 
col. In Sec. Ill CI we review some applications of the isomor- 
phism Known distance measures for quantum states can 
be used to provide distance measures for (trace preserving) 
CPM, which we will use in the following. Finally we extend 
the Jamiolkowski isomorphism in Sec. Ill DI to the multi-party 
setting and discuss some impUcations for the entanglement 
capabilities of CPMs. For sake of completeness the reader 
can find a review yii^ about the relation between the spec- 
tral decomposition of states and the Kraus representation for 
CPM in Appendix A and about the relation between the purifi- 
cation for quantum states and quantum operation in Appendix 
B. Note that the main properties of the isomorphism are stated 
in the form of short propositions with a consecutive number- 
ing (No. 1 - 13) , that is continued in the Appendix. 



A. The Isomorphism in the general setting 

Let and Ha' be two Hilbert spaces of finite dimen- 
sions dA — dimc(H^) and d^/ = dimc(H^/). With /Aa = 
A^(H^) and A4a' ~ A^(H^') we denote the corresponding 
matrix algebras over Ha and Ha' respectively, which con- 
tain the set of physical states (density matrices) Va C Ma 
and Va' C Ma' as (proper) convex subsets. Similarly, we 
will write A^(Hyi,H/i/) for the algebra of d'^ x dA matri- 
ces representing the linear maps from Ha to H^/ . Moreover 
let End {Ma -Ma') be the set of linear maps (endomor- 
phisms) between the algebras A4a and AiA' , which contain 
the physical operations CPM {Va — > Va' ) between the two 
quantum systems, i.e. completely positive maps (CPM), as a 
proper subset. In the following we will frequently consider a 
copy A of system A and use - after identifying and fixing a 
basis in and Ha' - the maximally entangled state 
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1 



i=l 



(1) 



on the composite system ig) Ha- 
The Jamiotkowski Isomorphism 

The map J : M (H^/ (g) H^) End {Ma ^ Ma'), that 
maps a matrix E of the matrix algebra over the composite 
system Ha' ® Ha to the linear map £ given by 

£{M) := d\ trAA [e^'^ M^] (2) 

(for any M G Ma — Ma) 1^ isomorphism ^ 0/, i-e. 
it is linear and bijective. For the inverse of J the matrix 
E <^ M {Ha' ® Ha), that corresponds to the linear map £, 
is given by 



E 



^ld^(Pt) 



(3) 



If the matrix E and map £ in correspondence are decom- 
posed with respect to the chosen basis |j) {j G N^^) on Ha 
(or H^) and |/) (i € N^^,) on Ha 



E 



A' A 



E E,m\')^'ik\(g>\j)^{l\ (4) 



£{M) = ^^m{m\i) N>-^'< 



(5) 



''a 



the Jamiolkowski isomorphism simply is 



£, 



k\ji 



dAE, 



ij\kl 



(6) 



From this relation between the coefficients the bijectivity im- 
mediately follows from the linearity of J' together with the 
fact, that M {Ha' <^ Ha) and End {Ma ^ Ma') both are 
linear spaces of dimension dA x dA' ■ The above result there- 
fore can be shown by deriving relation (|6} separately from 
Eq. (13 and from Eq. Q using the fact that 



A\t 



(7) 



holds for any M e Ma- 

The isomorphism also turns out to be an isometry 0]; If 
one uses the inner products 



{E,F) 



tvE^F 



(8) 
(9) 



for E, F E M {Ha' <E) Ha) and the corresponding maps 
£,T E End {Ma Ma'). "^^ can readily be seen that 



{E,F) = {£,T) 



(10) 



holds. Note that the corresponding norms ||M|] on Ma 
{Ma') and on End (^M^ Ma') are Euclidean {li 
norms). 

Eq. (|6jl shows, that the isomorphy does not extend 
to the respective compositions in M (Ha' ® H^) and 
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End {A4a ^ -Ma')- For example, the composition E o T 
of to maps would correspond to a matrix composition law 
{EC)F)ij\u = Yl,m-n Errn\knF„ij\ni which differs from the 
usual matrix multiplication. Nevertheless it also provides 
A4 {a A' ® Ha) with a semi-group structure @|. 

The effect of matrix multiplication by local 01 matrices 
is given by the following formula: Given a map £ and its 
corresponding matrix E and the matrices Bi,Ci G Ma', 
^2, C2 G Ma, then the transformed matrix 

E' = Bf (g) E^ Cf (g) (1 1) 

corresponds to a map £' with 

£'{M) Bi £ {BlMCl) Ci . (12) 

B. The isomorphism for quantum states and quantum 
operations 

Under which conditions on the matrix E e 
M (H^' ® H^) does the linear map £ E 
End {Ma ^ -M A') correspond to a physical operation, 
i.e. is a (trace-preserving) CPM ? This can be answered by 
the following results 

1. £ is Hermiticity presen'ing^^, iff E is Hermitian. 

2. £ is positivity preserving |0], iff E' is Hermitian and for 
all separable states F G V (Ha' H^) tr (EE) > 
holds. 

3. £ is completely positive iff i? is positive. 

4. f is a trace-preserving CPM jsl, iff _B is positive and 
trA'E^'^ ^ -jjIa holds. 

These results imply that the JamioHcowski isomorphism can 
be restricted to J : V (Ua' ® Ha) ^ End {Va Va') 
yielding a correspondence between trace-preserving CPM 
and states on the composite system of A and A'. In this case 
the isomorphism can be given a natural interpretation in terms 
of a teleportation protocol (without classical communication): 
In order to obtain the state E corresponding to a CPM £ ac- 
cording to Eq. Q the CPM £ simply has to be applied at the 
system A of the composite system in the maximally entangled 
state |<i>) (see Fig.0. Conversely, given the state E, the CPM 




FIG. 1 : In order to obtain the state E the CPM £ is applied to system 
A of the joint system of A and A, which is prepared in the maximally 
entangled state P,^'^. 



£ can be evaluated for an arbitrary input state p according to 
Eq. (2} as follows (see Fig|3- Considering the composite sys- 
tem consisting of parties A' and A in the state E together with 
the input state p at system A, i.e. the total state E^ ^ ® p^, 
the joint system AA is measured in a Bell basis containing 
the maximally entangled state P^^. With probability -k- the 
desired output state £{p) is then obtained at system A' . 

A' A^- -A 

• ^ 

^ ^ Projection |" ^ 

FIG. 2: Given the state E on the composite system A' and A, the 
CPM £ is evaluated for an arbitrary input state p by taking p as an 
input at system A. Then the joint system AA is measured in a Bell 
basis containing the maximally entangled state P^^. With probabil- 
ity the desired output state £ip) is then obtained at system A' . 
"a 

According to Eq. il2\ any operation on E, which is sep- 
arable w.r.t. the partitioning {A' , A), i.e. 

E' = J2 Bf ® CfE{Bf' ® Cf)^ , (13) 
j 

translates to a probabilistic application of combined opera- 
tions before and after the CPM £: 

£'{M) ^3 ^ (CjMC*) b] . (14) 

j 

In particular, the application of local unitaries or measure- 
ments to E on party A {A'} corresponds to the application 
of local unitaries or measurements before [after] the CPM £. 
On the other hand, not all separable operations can be imple- 
mented by localquantum operations and classical communi- 
cation (LOCC) 121 ■ Since only the measurement results before 
the CPM £ can influence operations performed afterward, we 
have to restrict the separable operations on side A and A' even 
to be local quantum operations and one-way classical commu- 
nication (1-LOCC) from party A to party A' . The separable 
operations in question thus correspond to the state 

where 

(i) CjCj ~ 1, i.e. the quantum operation C{p) = 
J2j C^jPC^j on party A is a trace preserving CPM; 

(ii) C is bi-stochastic C(l) = 1 and hence the correspond- 
ing CPM C(p) = J2j C^JP^j before the application of 
£ is also a trace preserving CPM, i.e. J^j ^j^J ~ 1' 

(iii) for each measurement outcome j on party A [before the 
application of £] the corresponding operation Bj (p) = 
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^ ■ BijpBjj, that is performed on party A' according to 
the classical information sent by A, is a trace preserving 

CPM,i.e. E,4^.. =1- 

Condition (i) and (iii) specify the notion of a general 1-LOCC 
protocol, that we consider in the following. In many cases 
such as for local projective measurements or for probabilistic 
applications of local unitaries, property (ii) follows from (i), 
but in general (ii) provides an separate condition, which re- 
flects the fact that before £ not C but C with transposed Kraus 
operators is applied. To simpUfy notations we will therefore 
consider those 1-LOCC protocols, that satisfy all three con- 
ditions. The above discussion indicates the two directions, in 
which one can try to manipulate £ with the help of the corre- 
sponding state E: 

(A) If one really has the above teleportation protocol avail- 
able in practice, any (non-local) operation on E can 
be considered in order to manipulate the corresponding 
CPM. 

(B) If the isomorphism is only a helpful theoretical tool, 
then one should only consider 1-LOCC operations on 
E in order to manipulate a given CPM, since these op- 
erations can be implemented by a coordinated applica- 
tion of operations before and after the evaluation of the 
CPM. 

We emphasize that only one direction of the isomorphism pro- 
tocol can be implemented with unit probability of success. 
This implies that the case (A) corresponds to a probabilistic 
modification of the CPM £ whereas case (B) gives rise to a de- 
terministic manipulation protocol. Case (A) is also equivalent 
to all protocols, in which one does not only allow arbitrary lo- 
cal operations before and after the application of £ (and there- 
fore the use of independent ancillary systems to perform these 
operations) but also to make an (arbitrary) ancillary system 
available to store quantum information. This information is 
obtained during the operations before the CPM £ and later 
used in the operations performed after £ ifioll . 

C. Distance measures for quantum states and quantum 
operations 

In the remainder of this paper we derive standard forms £' 
for some noisy CPM £, that approximates some ideal opera- 
tion T. A reasonable requirement for such a standard form is 
that it is also a considerably good approximation to the ideal 
operation. In order to assess this requirement some kind of 
distance measure between quantum operation is needed. As a 
first application of the Jamiolkowski correspondence we thus 
review the derivation 0] of distance measures for quantum 
operations from those for quantum states. 

Concerning the isometry properties discussed in Sec. Ill Al 
note that the Euclidean norm does not provide a proper dis- 
tance measure d for quantum states, since it does not obey the 
contractivity property, that is 

di£ip),£{a)) <d{p,a) (15) 



for all states p, a and trace-preserving quantum operations £ 
fTlll . This property expresses the physical condition that no 
quantum process should allow to increase the distinguishabil- 
ity of two quantum states. In the literature (see e.g. fyJT^J]3\) 
there are mainly two metrics fl4ll considered that also obey the 
contractivity property, namely: 

• trace distance: di{p, cr) := ^l/c — cr|ti-, where |Af|tr '■— 
tr (yMH'i'j is the trace norm: 

• fidelity-based distances, that are monotonically de- 
creasing functions of the fidelity F{p, a) := 

^^(VVP'^Vp) such as d2(p,cr) := ^1 - F{p, a); 
note that ^(1-0), cr) — (■0|cr|V') if p = |V'')(V'| is pure. 

By again using the Jamiolkowski isomorphism both distance 
measures for quantum states have a natural counterpart as a 
distance measure for quantum operation. Given two trace- 
preserving quantum operations £ and !F one can define the 
distance A{£,!F) = d{E,F) as the distance d between the 
corresponding states E and F, which is easily shown to yield 
a metric A on the set of trace-preserving quantum operations 
as long as d is a metric on the corresponding set of quantum 
states. Choosing d ~ di or d ~ d2 the corresponding distance 
measures Ai and A2 also have the following two properties, 
which seem to be reasonable requirements for any distance 
measure for quantum operations |5|| : 

• Stability (H: A(Id (g) £,ld® J") = A{£,J'), i.e. 
the distance measure of two quantum processes should 
not depend on whether they are considered to occur in 
an environment together with some unrelated ancillary 
quantum system; 

• Chaining Q: A(£i o £2,^1 ° ^2) < A{£i,Ti) + 
A(£2, -^2), i-e. for a composed process, the total error 
will be less than the sum of the errors in each individual 
step. 

Note that the evaluation of the above distance measures in 
practice requires some quantum process tomography. More- 
over both measures can be shown to have some physical inter- 
pretation in the sense of a bound to the average-case-error in 
function computation and sampling computation jsj. But in 
the following we will consider a slightly different application. 
Unfortunately the natural approach for defining error mea- 
sures for quantum operations by averaging the distances be- 
tween the output states, i.e. A{£,!F) :~ J dip d{£ , T {ip)) 
so far could not be modified in such a way that it would also 
fulfill the stability property. Nevertheless, for the case that one 
operation = [/ is a unitary operation, the average fidelity 

F{£,U):= j d4> {4'\UH{\i>){ii:\)U\iP) (16) 

has at least a plausible interpretation in terms of the average 
'overlap' between the two outputs U]^!:) and al- 
though it does not even define a metric. It was shown in 
lll6ll that this average fidelity F{£, U) is linearly related to 
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the 'Jamiolkowski' fidelity F{£, U) = F{E, \i}u)) by: 

F{£,U)d+l 



F{£,U) 



1 



(17) 



where \'4)u) denotes the pure state corresponding to the unitary 
operation U . 

In the following we will be interested in standard forms for 
noisy operations £, which approximate some ideal operation, 
that will be either the identity Id or some unitary U . These 
standard forms £' are obtained by different protocols, which 
might introduce additional noise to the operation £. Apart 
from simplicity of the obtained standard form £' , it should 
only differ in the same order of magnitude from the ideal oper- 
ation U or Id as the original imperfect operation £. The above 
argument shows that the 'Jamiolkowski' fidelity F{£,IA) or 
F{£,U) can in both cases be used to measure this distance: 
On the one hand the fidelity is related to a decent distance 
measure for quantum operation by a monotonic decreasing 
function. On the other hand for our applications the fidelity 
has a physical interpretation in terms of the average error in 
approximating an ideal (unitary) quantum operation U. In the 
following we therefore try to provide standard forms £' of 
noisy operations £, that have either the same or a slightly de- 
creased fidelity F{£', U) with the ideal operation as the orig- 
inal one (F{£,U)). 



D. The Isomorphism in the multi-party setting 

The Jamiolkowski isomorphism has a natural extension to 
multi-party scenarios, which are of special interest in quantum 
information theory. For this let the system A = (Ai , . . . , An) 
consist of N parties, each representing Hilbert spaces H^* of 
different dimensions cIa,, such that = H^^ ® . . .® H^" 
and dA ~ YiiLi ■ In order to keep the argumentation sim- 
ple we consider only CPM £, whose input and output Hilbert 
spaces are of the same type, i.e. = H"*i (g) . . . (g) H"*« 
withH^> ~ H^'. 




FIG. 3: In order to obtain the state E the CPM £ is applied to the 
systems Ai of the joint system A — {Ai, . . . , An), which are (lo- 
cally) prepared in the maximally entangled states P^'^'. 

The main point in extending the isomorphism to the multi- 
party setting is to choose the maximally entangled state 



to be the tensor product of the respective maximally entangled 
states 



|a>) 



AiAi 



(18) 



fc=i 



between the subsystem Ai and its copy Ai at each individ- 
ual party i = 1, . . . , TV. The maximally entangled state uj is 
therefore 



with Pf'^' = 1$)^^^' 



(19) 



l^j. In this notation the isomor- 
phism will have exactly the same form as stated above with 
the only difference that the maximally entangled state $ used 
in both directions now also respects the partitioning A = 

{Ai,...,ANy. 



£{M) 
E 



di tr 



AA 



^A'A pAA ^^A 



£^®UA{Pi^) 



(20) 
(21) 



E 




Projection 



Input 



FIG. 4: Given the state E on the composite system A' — 
(A'l, . . . , A'n) and A = (Ai, . . . ,An), the CPM £ is evaluated 
for an arbitrary multipartite input state p by taking p as an input at 
system A = [Ai, . . . , An)- Then the joint systems AiAi are (lo- 
cally) measured in a Bell basis containing the maximally entangled 

state p_^'^'. With probability the desired output state £{p) is 

"a 

then obtained at system A' . 



For the interpretation in terms of a teleportation protocol 
the preparation of the maximally entangled state P^^ and the 
corresponding Bell measurement can be performed locally at 
each party separately, since the entanglement present in P$ 
is only with respect to the systems Ai and their copies Ai 
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but not with respect to the partitioning itself (see Fig. |3] and 
Fig. 1^. For the index notation it is convenient to take the 
same formula as in Eq. (|6} 

^-ik|jl = dA ^'ijlkl , (22) 

but to consider the indices i, j, k and 1 as multi-indices, e.g. 
i = (ii, . . . ,ijv) e Nd^j X ... X Nd^^. It turns out, that 
many of the entangling capabilities of the CPM £ are directly 
related to the entanglement properties of the corresponding 
state E (see iHITtIi^: 

5. £ is separable w.rt. parties Ak and Ai (and therefore 
not capable to create entanglement between them), iff 
E is separable w.r.t. parties Ak and Ai. In particular, 
we find that the CPM corresponding to the tensor prod- 
uct of states E ® F simply is the tensor product of the 
corresponding CPMs £ ® T. 

6. For the partial transposition w.r.t. party A^ we 
have: 

{£{M)f'''^^£'[M^^^) (23) 
withi;' = E'^'-^'k'-"" . 

In particular, £ is PPT preserving w.rt. party Ak O, 
iff E is PPr w.r.t. the joint transposition of A'f. and Ak- 

7. The CPM £ can simulate another CPM T under 
SLOCC figil . iff the corresponding positive operator E 
can be converted into F by means of SLOCC. 

8. Two CPMs £ and T are equivalent under local uni- 
taries (LU), iff the corresponding positive operators E 
and F are LU-equivalent w.rt. the finest partitioning 
{Ai,A[,...,AN,A'r^) . 

9. The CPM £ can generate a state p of the composite sys- 
tem H"^ (Xi with non-zero probability of success, iff 
the corresponding positive operator E can be converted 
into p by means of SLOCC. 

Since the classification of pure states in bipartite and three- 
qubit systems under SLOCC is known in detail, the results 
can be transferred to the corresponding maps via No. 7 and 9 
(see iHtIi '). For further applications to purification, storage, 
compression, tomography and probabilistic implementation 
of non-local operations and its use in quantum computation 
we refer the reader to Ref . ll20ll and 1 2 1 1 . 

In the following we will mainly consider two-qubit gates, 
which are of special interest in quantum computation and 
quantum information. Note that in Ref. 0] and 1I20I1 (see 
Sec. Ill) it was shown for the case of two-qubit unitary opera- 
tions, how to modify the teleportation protocol to implement 
an arbitrary two-qubit unitary with unit probability of success. 
We will illustrate the results No. 7, No. 8 and No. 9 for these 
gates, namely for the 

. CNOT-gate: U^^ = |0)^(0| + \1)^{1\<J^ 

I^cnot) = (lOO)^^'l^o)^^' + 111)-^^' IVi)''''') (24) 



• Phase-gate: 17^^ (a) = e'^^^y^^y 

\il){oi)) = cos(q)|7/)o)'*'^ \ipo)^^ - is\n{a)\ip2)^^ l^a)"^^ (25) 

• SWAP-gate: U"^^ = |00)'^'^(00| + |01)'^'^(10| + 
|10)^^(01| + |11)'4^(11| 

Inswap) = |V'o)^^'|V'o)^^' (26) 

Note that I^/iswap) is a product state w.r.t. the partitioning 
AB' versus A' B but not w.r.t. the partitioning AA' versus 
BB' . In fact jV'swAp) has a Schmidt decomposition into 4 
(AA', i3i3')-product terms, whereas IV'cnot) and can 
be decomposed into 2 {AA' , SB') -product states. From basic 
facts ll22il about bi-partite entanglement for pure states it fol- 
lows from No. 7 that the SWAP-gate can simulate the phase 
gate and the CNOT-gate by means of SLOCC operation to be 
performed before and after the SWAP operation (but not vice 
versa). Moreover the CNOT-gate and the phase gate can sim- 
ulate each other under SLOCC for arbitrary a sjO, 2ti\. For 
the case of a = ^ they actually coincide up to some local 
unitaries iHtIEsII : 

U^B^ = (S) Ui C/^^ ( J) V^^ ® Vi with (27) 

This corresponds to the fact that the corresponding states 
IV'cnot) and are LU-equivalent w.r.t. to the partition- 

ing [A, A', B, B') (see No. 8), i.e. 

I^cnot) = Uf ® Uf ® (Fi^)^ ® (y/)^ m^)) . (28) 

According to No. 9 the SWAP-gate will moreover be capa- 
ble to create more entanglement than the CNOT-gate and the 
phase-gate, which can - up to SLOCC - create the same type 
of entanglement. 

III. STANDARD FORM FOR DECOHERENCE IN THE 
SINGLE-PARTY-SETTING 

In this section we apply the Jamiolkowski isomorphism in 
order to derive a standard form for an arbitrary decoherence 
process, that is described by some CPM £. We show that this 
standard form can be achieved by randomly choosing appro- 
priate unitaries to be performed before and after the actual 
CPM occurs. We first consider the case of a qubit system and 
then discuss a generalization to d-level systems. 

Let denote ctq = 1, ai = (t^, a2 = and 173 = the 
Pauh matrices. Note that starting with the maximally entan- 
gled state 1$)^^ = -i= (^|00)^^ + 111)^-^) we obtain a com- 
plete Bell basis (|V'o), IV^i), 1-02), IV'3)) simply by applying cr^ 
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locaUy on system A i.e. IV',)^^ = afi^M^ (^It/'o)'^'^ j . Thus 

any decomposition of a state E = J^ij in terms 

of the Bell basis corresponds to a canonical representation of 
the CPM £ in terms of Pauli matrices 

3 

£{p) = ^ E.jOipaj , (29) 

where the conditions on the matrix E = [Eij ) can directly be 
read off from the isomorphism, i.e. E must be density matrix. 
A case of particular interest in quantum information theory, 
especially in the study of fault-tolerance of quantum compu- 
tation, is when is a diagonal matrix. In many applications 
the corresponding CPM, the so called Pauli channel 

3 3 

S{p) = ^ E,a,pa, with E, = 1) , (30) 

describes some underlying noise model or decoherence pro- 
cess. This class contains for Eq ~ ^^^^ and Ei = E2 = 
E3 ~ the depolarizing channel (white noise) £p = 

pp +{l-p)^l, for Ea = i±£, Si = £^2 = and ^3 = ^ 
the dephasing channel £ p = pp + [p + a^pa^) and for 
Eq = E2 ^ E3 =: and El = the bit-flip channel 

£p=PP+ ^ {p + (TxPC^x)- 

We show now that the decoherence process specified by an 
arbitrary CPM £ as in Eq. (I29> can be transformed into a Pauli 
channel £' (see Eq. ( I30» with the same diagonal elements 
Ei = Eii. This can be achieved by a probabilistic but corre- 
lated application of one of the four Pauli matrices ct; before 
and after the actual noise occurs; 

1 ^ 

1=0 

In other words, by randomly choosing one of the four Pauli 
matrices with probability j to apply to a system before and 
after the noise process affects the system (e.g. some memory 
device) and ignoring the information about which Pauli matrix 
has been applied, an experimenter will actually (only) have to 
deal with noise of the form of a Pauli channel. The fact that 
the CPM £ can be brought to this form £' follows from the 
Jamiolkowski isomorphism used as in case (B) and the fact, 
that the corresponding state E can be diagonalized to E' by a 
mixing procedure, in which each of the local Pauli operators 
(jf (S> erf is applied with probability j. 

j^,AA' ^^^Y.at® of E^^ at ® af . (32) 

1=0 

The achieved standard form in Eq. ( I30l l can be further depo- 
larized, by considering the following three Clifford operations 
Qk = e*T'^'= = ^Jiak with k — 1,2, 3. Starting with a state of 



standard form Eq. (I30t one can in fact compute that 

1 ^ 

-Y,Qk® QlEQl ® Ql (33) 
fe=i 

= Ea\'4ja){'4ja\ + {Ei+E2+E3) x 

X {\A){^i\ + \^2){i'2\ + \i'3)m) ■ 

This means that by uniformly choosing one of the 12 unitaries 
Uki = Qk<^i for k = 1,2,3 and i = 0, 1, 2, 3 and applying 
Ul,- before and Uki after the application of an arbitrary CPM 
£ (see Eq. (|29j) the resulting CPM £': 

ki 

will be of the form of depolarizing channel 

f'(p)=a(/)p+(l-a(/))trpil (35) 

with / = E'oo and a{f) = '^^^^ ■ Note that a similar 
twirling procedure is also used in the recurrence protocol ll24ll 
for entanglement purification. Both depolarization procedures 
Eq. ( 13 U and Eq. ( I33> leave the state 1$) and thus the identity 
operation Id invariant. Hence the Jamiolkowski fidelity re- 
mains the same, i.e. F{£,\d) = _Eoo = F{£' Since the 
Jamiolkowski fidelity represents the noise level of the respec- 
tive operations £ and £' , both standard forms can be achieved 
without introducing additional noise to the system. 

Let us now turn to the case of general qudit systems with 
d = dimc(H"^) = dimc(H'^ ). Here the following complete 
basis of maximally entangled states can be chosen: 

■■= J2 + 0^|m>^' , (36) 

^ m—O 

where addition m + I and multiphcation fc • 7ti is meant mod- 
ulo d. Note that the Bell basis can be generated by acting on 
only one of the systems by means of unitaries Uki (generalized 
Pauli group) 

f/fc,|m) :=e^^'=-"|m + (37) 
out of the maximally entangled state |'0oo) = e.g. 

li^ki)^^' =U^i®ldA'[\%o)''^') ■ (38) 

Similar to Eq. ( I29l l the canonical form for an arbitrary CPM 
in terms of the generalized Pauli operators is 

£{p) = ^t^i^k-i' UkipUii, . (39) 

With respect to this Bell basis the corresponding state has the 
decomposition 

E= ^kLk'i' IMi^k'i'l ■ (40) 

kl,k'l' 
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By generalizing the depolarization procedure in Eq. ( 13 1> 
we can again diagonalize the state E and thus bring the corre- 
sponding CPM £ to the form of a (generalized) Pauli channel. 

Standard form: Pauli Channel 

By uniformly choosing one of the dP Pauli operators Uki and 
applying before and Uki after the application of an arbi- 
trary CPM £ (see Eq. |OP}j the resulting CPM £' : 

E Um£{uIpUu)uI (41) 

will be of the form 

£'{p)^Y.^kiUkipUkit (42) 

kl 

with E'f^i = EkiM- 

A proof of this statement can be found in the Appendix 
C. Whereas an arbitrary hermitian matrix E is described by 
— 1) real parameters, which in addition have to fulfill 
the constraints No. 4 (Sec. IIIBI i in order to correspond to a 
trace preserving CPM £, an arbitrary Pauli channels f can be 
described by only — 1 positive parameters E^i i25il . 

The number of parameters can even be decreased by con- 
sidering for the symmetrization procedure not only the Pauli 
group S := {±1, ±i} x {cro, ci, o'2, 0-3} but the larger group 

S' = {U^ ®U*^' \U €\]{d)} (43) 

of all local unitaries of the form ®U*^ . Since the group 
S' contains S, it has at most a smaller commutant. Whereas 
the commutant of S is the set of all Bell diagonal states, the 
commutant of S' is indeed t26ll only generated by the (orthog- 
onal) states P$ and 

(k, 1)^(0,0) 

In other words the set of states, that is invariant under 5', is 
determined by a fewer set of parameters. In fact, the states 
E' — ^{E) obtained by this 'twirling' operation 

V{E) = j [U ®U*)E{U^ ®U^)dU , (45) 

where dU denotes the uniform probability distribution on the 
unitary group U(d) proportional to the Haar measure, is de- 
termined by a single real parameter : 

E' = ./P* + (!-/) 7 (46) 
= a(/)P$ + (l-a(/))^l, (47) 

witha(/) = ^^fji. Note that the fidelity / = ($|£;'|$) = 
(^IPI^) (0 < / < 1) is left unchanged under the twirUng 
procedure T), since T) simply is a projection onto the subset of 
states invariant under this 'isotropic symmetry': 

V{E) = \x{P^E) + [d^ - 1) tr(7£:) 7 . (48) 



We remark that by partial transposition these isotropic states 
E' are in one-to-one correspondence I27il to the set of Werner 
states 101 ■ Since E' is a mixture of the maximally entangled 
state Pj, and a maximally mixed state, we find that the nor- 
mal form of the corresponding CPM £' is the (generalized) 
depolarizing channel 12811 : 

Standard form: Depolarizing Channel 
By uniformly choosing a unitary U £ U((i) and applying 
before and U after the application of an arbitrary CPM £ ( see 
Eq. (Oj the resulting CPM £' : 

£'{p) = j U£{U^pU)uUU (49) 

will be of the form 

(fc, 1)9^(0,0) 

= Q(/)p+(l-a(/))trpil (50) 

a 

with f ^ Eoo and a{f ) = . 

Since an isotropic state E' as in Eq. ( I47t is separable iff 
/ < i |26], we note that, according to No. 13. in Appendix 
A, the corresponding depolarizing channel becomes entangle- 
ment breaking at this point. 

Let us briefly address the question of possible practical im- 
plementations of the twirhng protocol described above. As it 
is shown in the Appendix C it is actually sufficient for the de- 
polarization protocol to uniformly choose some unitaries from 
a finite set of Clifford unitaries. 

To summarize we have shown that both standard fonns £' , 
the Pauli channel and the depolarizing channel, can be ob- 
tained by a random application of quantum operations ap- 
plied before and after the actual CPM £. These operations 
are chosen uniformly at random from a finite set of unitaries. 
Moreover we have seen that these depolarization protocols do 
not introduce additional noise to the system. 

IV. STANDARD FORMS FOR CPM IN THE MULTI-PARTY 
SETTING 

In this section we will continue the discussion of standard 
forms for noisy quantum operations. We will consider a par- 
titioning of the system A = (Ai, . . . ^Am) into N parties, 
which might be located at distant places. Any depolarization 
protocol that brings a given (non-local) CPM into its standard 
form therefore should be local w.rt. this partitioning. In the 
following we will consider an ideal operation I, that can only 
be realized imperfectly as a CPM £. We are now interested in 
the possible normal forms £' , into which we can transform £ 
by means of LOCC operation (w.rt. to the given partitioning), 
that are carried out before and after £ ll2^ is actually applied. 
If one is interested in the standard form for a map describing 
a given decoherence process itself, the ideal operation is the 
identity I = Id. Apart from the identity we will in the multi- 
party setting also consider the case, where the ideal operation 
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is some 2-qubit unitary operation I = U (U{p) = UpW), 
which can only be realized in form of some noisy quantum 
operation £. In contrast to the case discussed in the previous 
section, the locality requirements now impose rather severe 
constraints on the allowed operations to manipulate a given 
CPM. Note that the state / corresponding to the ideal opera- 
tion (identity or unitary) is pure. That is for X = U [T = Id] 
we have / = IV"?/) (V'al [I = respectively. Before we 

give an outline of this section we mention several aspects of 
the problem of finding such a standard form. 

• One can distinguish the two cases where only deter- 
ministic or also probabilistic transformations are con- 
sidered, i.e. whether it is possible to transform £ into 
the respective normal form £' in all of the possible mea- 
surement branches of the LOCC protocol or in at least 
one. 

• Closely related to this distinction is the question 
whether one uses the teleportation protocol directly as 
in (A) or indirectly as in (B), since a direct use of 
the isomorphism protocol in general has only a certain 
probability of success. 

• Firstly, one would like the transformation protocol T) 
(on state level) to leave the ideal operation invariant, 
i.e. V{I) = /. In this case the fidelity = \xIE' 
of the ideal operation with the transformed noisy oper- 
ation £' will be the same as the fidelity F{£,T) = trIE 
of the ideal operation with the noisy operation £. Since 
the Jamiolkowski fidelity with the ideal operation can 
be regarded as some kind of distance measure, the trans- 
formation will keep £ as close to the ideal operation as 
before. For the case of the ideal operation being the 
identity, the protocol V simply should be unital. On 
the other hand one might as well be allowed to sacrifice 
some fidelity with the ideal operation in order derive 
simpler standard forms. 

• The transformation protocol might bring any CPM onto 
its respective standard form {universal protocol) or it 
might be designed to transform a specific CPM into 
standard form . 

Note that most of the differences in these versions of the prob- 
lem only become important in the multi-party setting. This is 
mainly due to the fact that the depolarizing channel already 
provides a standard form for an arbitrary noise process, which 
can be achieved deterministically by a unital transformation 
and which is already specified by a single noise parameter /. 

We generalize the results of Sec.|III]in Sec. II V Al and de- 
rive standard form for decoherence processes (i.e. for the case 
T ~ Id) under the constraint that the underlying control opera- 
tion have to be local w.rt. the given partitioning. In Sec. llV Bl 
we discuss the case where the ideal operation is one of the uni- 
tary gates SWAP, CNOT or a phase gate with some arbitrary 
angle. In this section we restrict first to those depolarization 
procedures that are universal, deterministic and leave these 
unitary gates invariant. In Sec. llV Cl we also discuss the case 
where the fidelity / of the operation is decreased by a certain 



amount (i.e. additional noise is introduced) in order to obtain 
a simpler standard form describing the noise process, that is to 
reduce the number of required parameters. For gates locally 
equivalent to SWAP and CNOT, this leads to noise processes 
described by global white noise. A similar result is obtained 
for all phase gates, provided that one has control over switch- 
ing the noisy operation on and off at will. The noise is - in 
the worst case - increased by an order of magnitude. Finally 
we briefly mention the problems in Sec. llV d1 that occur when 
trying to transfer the techniques developed in Sec. IIV Bl and 
Sec. IIV Ci to the more general case of an arbitrary unitary op- 
eration as the ideal operation. 

A. Standard forms for decoherence in the multi-party setting 

1. Depolarization without sacrificing 

Let us now consider possible standard forms for noise oper- 
ations (i.e. ideal operation is the identity) in the multi-partite 
setting. Note that the twirling operation 2? used in Sec.lIIIIcor- 
responds to a projection into the space of states E over H'^ 
H-^' (recall that we chose dA = dimc(H^) = dimc(H^') = 
Y[f=i '^Ai )■ Thus it is straightforward to derive the corre- 
sponding standard form for the multi-party case, which is ob- 
tained after sequential application of the twirling operation 
T) locally at each party. Since the invariant group in ques- 
tion are (^.^j and {^^^ S'j^, respectively, the commu- 
tants of these groups (i.e. the subspace of invariant matrices 

in M {^^' ® H-^^ , onto which the projection <S)'t=i T^Ai 
projects) are just given by the tensor products of the commu- 
tants (subspaces) for each party. Hence, a probabilistic appli- 

cation of local unitaries gki = Qkih ^ ® ■ ■ ■® Qk^iiJ^ v^'\\h 
probability ^ = ■ ■ ■ -jj— diagonalizes E and gives the 
generalized multi-partite Pauli channel 

£'{p)= J2 E^iUkipUki, (51) 

which is again specified by c?^ — 1 positive parameters (Nd^ : ~ 
{0, . . . ,di — 1}). In the case of equal dimensions d the channel 
is determined by dA = d'^^ — 1 parameters. For a noise oper- 
ation on two qubits, for example, the corresponding standard 
form is given by 

3 

£' (p) = ® '^f P'^f''^ ■ ^52) 

As in the case of a single system, further depolarization is 
possible and hence a simpler standard form can be achieved. 
To this aims one performs a complete twirl over the larger 
group S'. The result of this twirl is that one projects E into 
the set of states of the form 

E ^k7k, (53) 

keN2X---xN2 
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where 7k = 7^^^ ( 



■'^^k" each party Ai 7^' denotes 

P$ or 7i =7 spanning 



one of the two orthogonal states 70 
the respective isotropic subspace. If E' is decomposed w.r.t. 
the (non-orthogonal) basis (70,71) = (P*, 32— 1), the cor- 

responding CPM £' has a natural interpretation in terms of 
different white noise factors. To be more precise let us con- 
sider the example of two qudits with d = (Ia^ = ^^2- In this 
case the state E' is of the form 

E' = aoo Pi' ® Pi' +aoi Pi' A, 

^ .P^= +aii jjlAi ® ^1^2 (54) 



Since {5q, 61) = (P$ — 7, d?^) is a dual basis for ^1), 
we have that 



"00 


= tr 


Ml 
"0 










Eqi Eio 




d2 - 1 d^- 


- 1 


aoi 


= tr 


Ml 
"0 








d^Eoi 


d^En 








- 1 


(d2 - 1)2 




Qio 


= tr 


'st' 








d^Ew 


d^En 








- 1 


(d2 _ 1)2 




ail 


= tr 


\5t' 







(d2 - 1)2 



;En ■ (55) 



(d2 - 1)2 

For the corresponding normal form £' of the CPM £ we obtain 



£'{p) = aooP + "01 PAi (X) 3IA2 + "10 3IA1 •X' PA2 

d d 



1 , , 
+"11 ^tr(p)lAiA2 



(56) 



where pA^ — trA2{p) and PA2 — ^Ai{p)- The second and 
third term correspond to white noise introduced locally at 
each party, whereas the last summand introduces global white 
noise. Note that some of the coefficients ago, "01 or can 
be negative. In the case of equal dimensions at each party 
{dAi = d) we have reduced the number of parameters from 
initially Oid^'^) for a general CPM £ over ^2^ _ 1 for a 
general multi-party Pauli channel to 2^ — 1 parameters (in- 
dependent of the dimension d) to describe different types of 
multi-party white noise. 

So far we have only considered twirling protocols, that were 
deterministic, made direct use of the isomorphism and left the 
ideal operation invariant, namely the identity operation. But 
can we further reduce the number of parameters for a differ- 
ent standard form, which is achieved by a LOCC protocol, 
that has only a certain probability of success or that makes in- 
direct use of the isomorphism or allows to sacrifice some of 
the initial fidelity E'oo with the ideal operation ? To be more 
precise, is it possible by weakening one of these conditions to 



achieve e.g. just a global white noise channel 

£'{p)^PP+il~p)^l, 



(57) 



i.e. the case, for which also all the diagonal elements vanish 
except aoo and an? 

If one allows for a direct use of the isomorphism protocol 
and thus to perform local Bell measurements in the basis 
this is certainly possible, since the ideal operation Id is local 
w.r.t. any partitioning and the coefficients E'ki in Eq. ( 15 1> can 
deliberately be adapted without increasing the noise level by 
simply 'twirling' all components _Eki but E'oo into ^1. Note 
that this procedure gives rise to a probabilistic LOCC proto- 
col, since a direct use of the Jamiolkowski isomorphism can in 
general not be achieved with unit probability of success (see 
Sec. Ill B> . As discussed in Sec. Ill Bl this use of the isomor- 
phism will therefore in general not be of great interest for all 
applications, in which one would like to work with a standard 
form of a given CPM rather than with the CPM itself. 

Using the isomorphism only indirectly, we have to restrict 
the transformations to the class of SLOCC operation that are 
also local w.r.t. to {Ai,Ai). Note that one cannot increase 
the fidelity / with the ideal operation X = Id by means of 
any physical protocol, whenever / corresponds to the largest 
eigenvalue in E Dili, e.g. Eq. dsTt with Eqq > i?ki (^ k = 



1 = 0). In all these cases any SLOCC operation Ci (g) C2 , 
that is contained in the transformation protocol and that does 
not leave |$) invariant, will cause a decrease in the fidelity 
/' < / of the respective standard form £'. Thus any uni- 
versal transformation protocol, that yields a respective stan- 
dard form for all CPM without decreasing the fidelity, con- 
sists in a probabilistic application of operation Cf (g) with 



Cf (E) Ci\^) = 1$), i.e. Ci®C2=C® {C-'^V for some 
invertible matrix C I32I1 . Since all such transformation will 
also leave all 7k in Eq. J53l l invariant, the respective standard 
form cannot contain fewer terms than the multi-party white 
noise channel. In this sense the above twirling procedure al- 
ready yields a standard form, that is optimal among all forms 
achieved by some universal protocol, that only use the iso- 
morphism indirectly and does not sacrifice any fidelity with 
the ideal operation. 



2. Depolarization by means of sacrificing 

In the remainder of this subsection we will show how to 
design twirling protocols that bring a specific CPM into the 
standard form of global white noise by introducing additional 
noise. The procedure described below does therefore satisfy 
neither the universality property nor the no-sacrificing condi- 
tion, but it will - especially for the many party case - signifi- 
cantly reduce the number of parameters of the standard form 
to a single one. By applying the above universal depolariza- 
tion procedures, we can start with considering only states E, 
that are already in isotropic form (see Eq. ( I53t ). 

We will illustrate the procedure for the case of two qubits. 
Generalization to the multipartite case and higher dimensions 



12 



are straightforward. For such two-qubit maps, the corre- 
sponding state E (after applying the universal depolarization 
protocol described in the previous section) is given by 

E = eooP^:pI^ + E + 



i=i 



E + E pf: ' (58) 

2—1 ijj — 1 



where = !?/'«) ("04 1 is the projector onto one of the Bell 
states \'4)i). In the following we will collect the parameters 
in a vector E = (E'oo, £'oi, -Bio, The fact, that E 

corresponds to a trace preserving CPM then simply reads (i) 
E > component-wise (i.e. Eij > for i, j ~ 0, 1) and 
(ii) Af(E) Eao + 3 {Eqx + Eio + SEn) = 1. We now 
consider the following type of depolarization 



(59) 



'3 'Pio. 

3=1 i=l 



-Pii E ^ '^t'E'^f' ® ^. 



In a similar notation as before V corresponds to a trace pre- 
serving CPM iff A^(p) = 1 and p > 0. It is straight for- 
ward to calculate, that the resulting state E' = ^{E) is again 
in isotropic form Eq. ( I58l l with new parameters E'^j, that are 
given by the linear transformation 



E' = D[E].p, 
where the matrix D[E] depends on the initial state E: 

(Eqq 3Eqi SEiQ 9Eii 

Eqi Eqq 2Ecii 3Eii 3{EiQ -\- 2Eii) 

EiQ 3Eii Eqq -\- 2EiQ 3{Eqi -\- 2Eii) 

Ell EiQ + 2Eii Eqi + 2Eii Eqq + 2{Eqi + £^io + 2Eii 



(61) 



Our goal is to bring E into a form, that corresponds only to 
global white noise (see Eq. (|57}), i.e. E'^^ = E[q = \E[^. 
Since we require E' > and 7V(E') = 1, E' should be of 

the form E'(/') = [f', with the new fi- 

delity /' G [0, 1]. Thus we look for solutions p to the lin- 
ear system Eq. ( I60> for the specific choice of E'(/'), that 
additionally fulfills the constraints A^(p) = 1 and p > 0. 
One computes that A^(E') = iV(p)A^(E) and we can there- 
fore omit the trace preservation condition A^(p) = 1, since 
we already require (chose) E and E' to be trace preserv- 
ing. Using A^(E) — 1 one can compute for the determi- 
nant det(D[E]) = rst, where r = 1 - 4:{Eoi + 3Eii) , 
s = l-4:{Ew+3Eii) andt ^ l-4(£'oi+-Bio+2-Bii). Thus 
the linear system Eq. ( I60> will definitely have a unique solu- 
tion, whenever max (£'oi, i?io, < (or alternatively 
the initial fidelity / = Eqo > 4max (i^oi, i^io, £^ii))- Let 
us consider a fixed vector E with det(D[E]) ^ first. For 
the corresponding CPM £ we want to design a standard form 
£' with maximal fidelity /'. In contrast to the standard forms 



discussed so far the depolarization process V, which trans- 
lates into applying Pauli operators before and after the actual 
CPM £ occurs, will be specifically designed for the given ini- 
tial CPM £, since the corresponding probabilities are given 
by the unique solution p(/') = D[E]"i • E'(/'). Note that 
each of the inequalities pij > is hnear in /', i.e. of the form 
flij /' + bij > 0, where the coefficients are 



1 

" 6 
1 

18 
1 

18 
1 

54 



o-oo 



a-oi 



aw 



1 1 4 
- + - + - 

rst 

3 1 4 
s r t 

3 1 4 
rst 

4 3 3 
t r s 



boo 
boi 
bio 
bii 



1 

48 

1 
144 

1 
144 

1 
432 



1 5 

1^ t 
5 1 

'l-r 
5 1 



3 3 5 

27 

rst 



Depending on the signs of 7^ the constraints are thus rep- 



resented by intervals starting or ending at /j' 



the corresponding condition is either always or never satis- 
fied). In the following we will discuss the (complete) positiv- 
ity condition p > in terms of the parameter (r, s, t) instead 
of {Eqi, Eio, Eu) since they are linearly related. Because of 

< Eij < 1 we generally have —15 < r,s,t < 1. The 
back transformation is given by Eqi ~ (1 + 3s — r — 3t), 
Eio = (1 + 3r - s - 3t) and = -^{1 + t-r- s), 
which implies £^00 = (1 + 3(r + s + 3t)). 

Let us consider the situation, in which the initial CPM 
£ is close enough to the ideal operation, such that 
max {Eoi, Eio, Ell) < Jq- Then there exists a unique so- 
(60) lution with 1 > r,s,t > 0. Moreover a further restriction 
of r,s,t to the interval will ensure aoi,aio,aii < 

and /qq < /oi,/io>/n- Thus for a fixed initial vector 
E we are left with the three constraints /' < /oi,/ioi/ii 
(and /' < 1) and the maximal achievable fidelity is /^^^ = 
min (1, flfi, f[„, f[i). Although restricting Eio, Eqi, Eu to 
the interval [0, ^] will be sufficient to guarantee that r,s,t E 
]|,1], not all r,s,t in the interval correspond to 

-£'01 J Eoi, Eoi G [0, 1]. A minimization of /^^^ for r,s,t £ 

1 1 , 1] therefore yields only an upper bound to the fidelity de- 
crease /i^^x (versus /) or the increase of 1 — /^ax (versus 1 — /), 
which represents the noise level of the CPM. A numerical 

minimization in this region shows that the relative fideUty 

is at least 37.14% and the relative noise level is at most 

increased by a factor 5.5. 

Note that for an initial fidelity / > y| we have < 
Eio, Eoi , Ell < 3^- In other words any decoherence process 
on two qubits, that introduces only little noise ( i.e. / > y|), 
can be brought into the form of global white noise by increas- 
ing the noise level by a factor less than 5.5. This standard 
form is achieved by application of local Pauli operations, that 
are chosen randomly according to some probability distribu- 
tion specified by the parameters py (/4ax)- Thus the protocol 
is specifically designed for the initial form of the decoherence 
process (more precisely it depends on the vector E, that is ob- 
tained after the decoherence is brought into the form Eq. (15 8> 
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by the methods described above). 

If the initial CPM does not belong to the region with < 
Eio, Eqi, Ell < ^ similar derivation can be applied. The 
constraints a^/' + bjj > again determine, whether a stan- 
dard form can be obtained in this way and how much fidelity 
has to be sacrificed in order to achieve the normal form with 
E'(/'). Moreover a generalization to the case of d-level sys- 
tems with d > 2 and to the multi-party setting with N > 2 
can be developed along the lines of the previous discussion. 
Note that for increasing N, although the achieved standard 
forms will also be global white noise and thus be specified 
by one parameter only, the derivation and the transformation 
protocol itself will become more involved, since the number 
of parameters Pi-^...iM will be 2^ and thus increase exponen- 
tially. 



B. Standard forms for noisy unitary operations 

We now turn to standard forms of noisy operations where 
the ideal operation is given by some unitary operation U. 
We concentrate on two-particle operations and will illustrate 
our approach with help of several examples, including gates 
which are up to local unitary operations equivalent to the 
SWAP gate, the CNOT gate and a general phase gate with 
arbitrary phase a. We will show that one can depolarize these 
noisy gates to standard forms with a reduced number of pa- 
rameters, without changing the fidelity of the ideal operation. 
To this aim, we decomposes a CPM £ into a unitary part U 
and some remaining (orthogonal) part £' (where £ is in gen- 
eral no longer a CPM), i.e. 

£p = fUpU^ + (1 - f)£p, (62) 

and both, / and £, are determined by the Isomorphism. We 
have that 

/ = i'^ulEliJu) (63) 

E ^ ^-f\^-)(^-\, (64) 
1-/ 

where E is the operator corresponding to the map £ and / 
specifies the initial fidelity of the operation U. It is not neces- 
sary to make such a decomposition, however in this notation 
it is immediately evident that only the noise part, namely £ 
is altered by the depolarization procedure. We remark that 
tr(£') = 1, however E may have negative eigenvalues and is 
hence not a density operator Nevertheless, we can formally 
decompose E (and thus £) into these two parts. We will show 
that one can depolarize the map £ to 

£'p ^ fUpU^ + (1 - f)£'p, (65) 

where £' is a (generally non-positive) map of certain standard 
form, specified by a few parameters. Clearly, the total map £' 
remains completely positive. Note that the depolarization of 
£ takes place in such a way that the weight of the ideal oper- 
ation is not altered. In particular, if the operation is initially 



noiseless (i.e. / = 1), it will remain noiseless after the de- 
polarization. This is achieved by considering depolarization 
processes that leave the unitary operation U (or equivalently 
the state \ 'ifu) when considering the operator E corresponding 
to the operation £) invariant. The number of required parame- 
ters and the expUcit form of £ depends on the ideal operation 
U, as the group of local operations that leave U invariant is 
determined by the structure of the state I'I'fy). 

1. The noisy SWAP gate 

In this section we determine a standard form for noisy 
SWAP operations. The ideal d-level SWAP operation 
is defined via its action on product basis states, namely 
f/swAP 10^1.7)'' = Ij)'^^)^, where {|fc)}fe=o,i,...,d-i is a ba- 
sis of H = C*. The state i?swAP = Inswap) (^'swafI cor- 
responding to INSWAP is specified by (see Eq. i26i ) 

Inswap) = (66) 

Consider the mixed state E describing — via the 
isomorphism — a noisy SWAP gate. We have that all 
operations of the form J/^ (g) U*^' (g) 1/^ ig) V*^' leave 
Inswap) invariant and hence can be used to depolarize E. 
This implies that we can essentially use the same depolar- 
ization procedure as in the case where the ideal operation is 
given by the identity (see Sec. lIV Al . only the role of particles 
A' and B' is exchanged. This implies that the resulting 
standard form can again be interpreted as a local and global 
white noise processes with three independent parameters, that 
occur before the application of an ideal SWAP operation, i.e. 

£'{p) = ;7swAp2?(p)t/swAP with 

11 

Vip) = aoa p + aoi Pa® -;iB + aio -1a Pb 

a a 

1 , , 

+q;ii ^tr(p)lAB . (67) 

Note that the parameters a^; are again given by Eq. ( I55L 
where E^i are the coefficients in a decomposition Eq. ( 15 3t 
of E according to the basis 

7k 6 { P^^' ® Pi^' , Pi^' ® 7^-^' , 7^^' ® , 7-^^' ® 7^^' } 

with 7 = ^2^1 {1 ~ P<s>)- In particular by the twirling 
procedure the Jamiolkowski fideUty remains the same, i.e. 
F{£, ?7swAp) = Fi£' , t^swAp)- 

2. The noisy phase gate and CNOT gate 
In this section we consider the unitary operation 

f/^B(«)=e-^"<^<, (68) 

for arbitrary angles a. Up to the local unitary operations, 
U (a) is equivalent to a controlled phase gate, while for a = 
7r/4, i7 (a) is equivalent to the CNOT gate (see Eq. j27li). i.e. 

t^CNOT = ® C/f C/(7r/4)l// ® V^^. 
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We will obtain a standard form for noisy operations, given 
in the ideal case by U{a), by depolarizing the corresponding 
CPM £. The depolarization takes place by applying appropri- 
ate random local unitary operations that leave the state 

= cos(a)|?Ao>AK'o>B - isin(a)|?/>2>A|'02>B, (69) 

invariant (up to an irrelevant phase), where |^'q,)(^'q,| is the 
state corresponding to U{a) via the Isomorphism, and 

1^',) =l®aj|$) (70) 

are Bell states. Note that such a depolarization procedure for 
U{a) automatically provides a depolarization procedure for 
all operations that are local unitary equivalent to U {a), lead- 
ing to a standard form with the same number of parameters for 
these noisy gates. The depolarization procedure simply has to 
be adopted according to the local unitary operations. To be 
specific, consider for instance the noisy [/(7r/4) gate and the 
noisy CNOT gate. If W^-^ is a local unitary operation that 
keeps |^'7r/4) invariant , then the operation 

{uf)^<E){Ui')^(g)V*f(g)V*^ 

keeps I^'cnot) invariant. That is, one obtains a depolariza- 
tion procedure for the noisy CNOT gate from the depolariza- 
tion procedure for the C/(7r/4) gate by replacing each unitary 
operation W by W'. 

We now present an explicit depolarization procedure for the 
noisy U{a) gate, described by the CPM £, with arbitrary a. 
We will consider the depolarization of the corresponding state 
by means of 4-local operations. We remark that any sequence 
of depolarization steps can be translated into a single step 
with multiple possibilities by considering all possible com- 
binations. Such a single step procedure can then be translated 
to appropriate random operations applied to the system before 
and after the application of £ and hence to depolarize the cor- 
responding map. For notational convenience, we define the 
four-qubit states 

|vI/,,)^^'^^'^|V.,)^^'®|^,>^^'. (71) 
Given an arbitrary CPM £ specified by 

3 

£p= ^ KjMO-iO-jPo-kO-i, (72) 
the corresponding state E is given by 

3 

^= E (73) 

where Xij^ki ~ A^; We define two-qubit unitary operations 
UM,Vhy 

U = (iay) (X) {i(Jy), 

U = (Tx®(J^, (74) 



The action of these operations on Bell-basis states can 
be readily obtained and one finds that U,U introduce rela- 
tive phases between the Bell states, while V exchanges two of 
them. To be specific, we have , 

W{|Vo),|^/'l),iV'2),|^/'3)} = {|^/'o),-|Vl>,!^2),-|^/'3)}, 
W{|^Ao),|^Al),IV'2),|^/'3)} = {|^/'0),|^l),-|^/'2),-|^3)}, 
V{|Vo),|^Al),l^2),|^/'3)} = {\^o),-\^'.i)M2)Ml)}. (75) 

All local operations of the form Ia^b, ^aI^b, I^a^b, 
UaUb, UaUb, IaVb, VaIb, VaVb keep the state 
(and the fidelity / = {'^ a\E\^ a) of the ideal operation) in- 
variant and can thus be used for depolarization. 

We decompose £ into the unitary part U (a) and the remain- 
ing noise part £ (see Eq. (I62» and consider the corresponding 
(non-positive) operator E (see Eq. (|S3}) in the following, 

3 

^= (76) 

i,j,k.l=Q 

We randomly apply IaIb, '^aUb, Ua'^b or UaUb, each 
with probability 1 /4, which leads to an operator 

E' = -^(e + U^EiU^y + U'^EiU'^y + U^U'' E . 

One finds that E' is of block-diagonal form with coefficients 
X'^j^ku that fulfill X'^^j^i = whenever {i mod 2) ^ [k 
mod 2) or (j mod 2) ^ [l mod 2) and remain invariant 
otherwise. This follows from Eq. i75\ . as U introduces a phase 
(—1) for Bell states with {i mod 2 = 1) while states 
with even parity (i mod 2 = 0) remain invariant, which 
results in the cancellation of the corresponding off-diagonal 
elements. Thus only elements A^^ ^.^ with (i mod 2) = (fc 
mod 2) and (j mod 2) = (/ mod 2) remain, which can be 
grouped into four 4x4 blocks Tab with a = {i mod 2) = (fc 
mod 2),b ~ {j mod 2) = {I mod 2). For instance, Fqi = 

Ei,fc6{0,2}u,i6{l,3} Ay,fe(|*y)(*fci|- 

In the following, we consider the depolarization of the sub- 
spaces Fat separately. We start with Fqo, which is spanned 
by the states {|*oo), |*02), |*2o), |*22)}- By randomly ap- 
plying IaIb or UaI^b with probability 1/2, we find that 
the resulting operator Fqq has coefficients Aqq q2 — Aqq 20 = 
Aq2 22 ~ 0, while the other coefficients remain invariant, i.e. 

Aqo^OO = AoO.OO; Ao2^02 = Ao2,02; A20,20 ~ A20,20; 
A22,22 = A22,22; Aoo_22 = Aoo,22 ; Ag2,20 = Ao2,20- (77) 

We thus find that Foo is of the form 



Aoo,oo 








^^00,22 


\ 





^02,02 


\' 

-^02,20 








^02.20 


-^20,20 







V A'oo,22 








\l 

•^22,22 





(78) 



which are 8 independent real parameters as XijM = A^; jj. 

The effect of these (random) operations on the other sub- 
spaces Foi,Fio,Fii is similar, i.e. the corresponding off- 
diagonal term vanish. However, in these subspaces further 
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depolarization is possible. Consider Fqi which is spanned by 
the states {|*oi), 1*03), l*2i), |^'23)}- Applying randomly 
IaIb or IaVb with probabihty 1/2 leads to coefficients 



^01,01 

^^21, 21 



A' 



03,03 



(Aoi,01 + Ao3,03) 



^^23, 23 - 2^'^21^21 ^ A23.23) 



-An- 



1 



(A 



01,23 



.23 ' 'V^.zi 2 

This can readily be seen by using that 

1aVb{|*oi), 1*03), 1*21), 1*23)} 

= {-|^'03),|«'0l), -1*23), 1*21)}. 

Thus we find that Fqi is of the form 



r' — 

^ 01 — 



^03,21 ) 



We remark that the depolarization process described in this 
final step leaves the subspaces Too, Foi, Fio — which were al- 
ready depolarized earlier — invariant. The final depolarized 
CPM is is specified by (8 + 4 + 4 + 2 - 1) = 17 real param- 
eters (where the (-1) results from the normalization condition 
tr(i?) = 1) and is of Block-diagonal form. The coefficients 
(79) Kj ki given by Eqs. I77I79I83I84I and are zero otherwise. 
This leads to the standard form, 

£p = fU{a)pU{a)^ + (1 - /) E Ay, « a, a, pa, a,, (86) 



01.01 








Aoi.23 \ 





^01.01 


\' 

^01.23 








, ,* 

01,23 


^21,21 





I* 

01,23 








A21,21 / 



(80) 



(81) 



and is thus described by 4 independent, real parameters 
(Aoi,oi A2i,2i are real, Aqj 23 is complex). 

Similarly, by randomly applying 1a 1b or Va^b with 
probability 1/2, one depolarizes the subspace Fio — spanned 
by the states {l^-io), l^so), 1*12), 1*32)}— to the form 



F' 



10 



/ A'lo.io 
A 


\ A'io,32 



A'io.32 \ 

10.10 ^10.32 

^ 10,32 ^^12, 12 'J 

X[,^,, J 



(82) 



where 



A', 



10,10 — 30 — 2^^10,10 + A3o,3o) 



Ai2,12 - A32 32 - -(Ai2,12 + A32,32) 



(83) 



A', 



10,32 - ^A3o,i2 — 2(^10,32 - A3o,i2), 



which is again described by 4 independent, real parameters. 
Finally, the subspace Fn — spanned by the states 
1*13), 1*31), 1*33)} — can be further depolarized by 
randomly applying one of the operations l^ls, IaVb, 
Va^b oiVa^b with probability 1/4. One finds that 

^^11,11 — -^11,11 ~ ^^13, 13 ~ ^31,31 ^ '^33,33 ^ 
— 7(Aii,ii + Ai3,i3 + A3i,3i + A33,33) 



A'll, 



33 



-A' 



1 



13,31 ~^m>^ll.33)~'fiiXu,3l)) (84) 



where 5R(.t) denotes the real part of x. Thus Fn is described 
by 2 independent, real parameters and is of the form 



F' — 
-■^ 11 — 



11.11 








Ail.33 \ 





A'11,11 


\' 

^11,33 








^^11,33 


All, 11 





11,33 








Aii,ii / 



(85) 



where/' = (*„|£;|*a) = i.e. the fidelity of 

the ideal operation remains invariant. To summarize, we can 
achieve the following standard form: 
Standard form for the Phase Gate 

By uniformly choosing one of the unitaries Ukfrom Ui -1^2 -1^3, 
where 



(87) 



Ui = {IaIb, e-*^"« Is, lAe-'^^^y , e" 
U2 = {I^Ib, o-fo-f } 

^3 = {I^Ib, O-^Ib, l^O-f , <^y<^y} I 

and applying Ul before and Uk after the application of the 
noisy phase gate £ the resulting CPM £' is of the standard 
form £'{p) = EL-.fe.i=o KiM^i^iP^^kcri with 



E' 



/no 




V 





F' 



01 



F', 



where 



10 



(88) 
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with the following choice of basis 



B 



{ 1*00) 
l^oi), 
l^io), 
l*ii>, 



1^-02) 

|*03>, 
l*30>, 
1*13), 



1*20) 
l*21>, 
1*12), 
l*31>, 



1*22), 
l*23>, 
1*32), 
|*33> } 



(89) 



and the parameters a,d,b,b,c,c,d,d,e,e € K and 
u,v,w,x G C. This depolarization does not increase the 
noise level, i.e. /' = /. 



3. The CNOT-type gate 

For certain values of a, further depolarization is possi- 
ble. In particular, we consider a = 7r/4, i.e. the opera- 
tions U{'k/A) which is local unitary equivalent to the CNOT 
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gate. In this case, the state |^'7r/4) is a maximally entangled 
state (with respect to systems A, B), which remains invari- 
ant under a larger set of local unitary operations than a non- 
maximally entangled state \^a)- In particular, we consider 
the unitary operations 



W = 1 (g) {iay), 

VV = Oz ® (Jx- 

which act on Bell states as follows 



(90) 
(91) 



W{|7Ao),|V'l),|V'2),|V'3)} = {i\i>2)M3),i\^l^0).~\^l)}.. 
W{|^0),|^1>,|^2),|^3)} = {-i|^2),|^3>,i|^0),|Vl)}. 

The operation WaV^b leaves the state |^'t/4) — up to an ir- 
relevant global phase factor (— i) — invariant. Note that this is 
not true for Y^a) with a ^ n/A. We take the standard form 
Eq. j86t as initial map, and apply randomly either Ia^b or 
VVaWb. One finds that the resulting operator E" is signifi- 
cantly simplified and is described by 8 independent, real pa- 
rameters. We denote the coefficients of E" by Hij^u- 
To be specific, for Fqq we find 



M00,00 — Ai22,22 - 2 (-^00,00 + ^^22, 22); 

M02,02 = M20,20 = 2 (-^02,02 + ^^20,20); 
Moo, 22 = *3'(^00,22);M02,20 = *^(-^02,2o)' 



(92) 



where 3(a;) denotes the imaginary part of x, i.e. 
i'^{x) = {x ~ x*)/2 and we thus have 4 real parameters. 
This follows from Wa VVb{|*oo), |*02), |*2o), |^'22)} = 

{1*22), -|^'20), 1*02), -|*00)}. 

Similarly, for V'^^ we find 



MOl.Ol = M03,03 = /^21,21 = M23,23 
A^01,23 = -A^03,21 = 3f^(-^01,23)' 

while F'/q simplifies to 

MlO.lO = /^30,30 = Ail2,12 = M32,32 
A^10,32 = —^30, 12 = 3fi(Aio.32), 



2 (-^01,01 + '^21,21); 



(93) 



^(•^10,10 + ^12,12); 



(94) 



where we have 2 real parameters in each case. 
Finally, for F'/^ we have 



Mil, 11 

Mil, 33 



Ml3,13 — M31,31 
-Ml3,31 = 0, 



M33,33 



^11, 111 



(95) 



which is a single, real parameter. It follows that the standard 
form for the depolarized gate U{'k/A) is given by 

£"p = /[/( J)mC/( J)^ + (1 - /) E l^^,M<y^<yJP<yk<Jum 

ij.kl 

where the coefficients jlijM defined in Eqs. ( I92L ( I93l l. 
( I94L ( I95> and are zero otherwise. Note that the fidelity of the 
ideal operation L/(7r/4) remains invariant. 



Standard form for the CNOT-type Gate 

The total state E is thus of the form 
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or equivalently 
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(97) 



vu 





where we use the basis B in Eq. H891 and a = Moo, 00. b 
M02.02 etc. are all real parameters. 



C. Standard forms by means of sacrificing 

While the standard forms for the general U{a) gate or the 
CNOT-type gate U{tt/A) are already relatively simple (as 
the number of relevant parameters is significantly reduced, 
namely from 255 to 17 or 8 respectively), for many practical 
applications a further simplification might still be desirable. 
If, for instance, one would like to analyze error thresholds for 
processes involving several particles and/or operations (as is 
e.g. the case in fault tolerant quantum computation or entan- 
glement purification with imperfect means) where noisy oper- 
ations are described by these standard forms, the correspond- 
ing CPMs are still rather complex. 

In this section we will provide such a further simplifica- 
tion of the corresponding noise process, where we find that 
in many relevant cases a single parameter is sufficient and the 
noise process can be described by (global) white noise. In 
contrast to the previous depolarization procedure, here the ex- 
act form of the noise process (equivalently the corresponding 
state E) has to be known. This may e.g. achieved by per- 
forming a process tomography of the CPM resulting after the 
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universal depolarization protocol describe in the previous sec- 
tion (note that only the knowledge of the depolarized map is 
required). In addition, the fidelity of the ideal operation is no 
longer conserved but decreased by a certain amount. That is, 
by "sacrificing" a (small) amount of the fidelity of the opera- 
tion, one can modify the resulting noise process in such a way 
that one obtains a very simple standard form. This is done by 
transferring weight from the ideal operation to the noisy part 
in an appropriate way and hence tailor the noise process. 



1. The noisy SWAP gate 

Let us consider a noisy SWAP gate in the two qubit 
case. For a specific noisy SWAP operation with suffi- 
ciently large Jamiolkowski fidelity / = F{£, Us^^p^^) = 
(V'swapI-E'IV'swap) > yI a depolarization procedure can be de- 
signed that brings the noisy operation £ to the standard form 



£'(p) = /'C/swAPpf/lwAP+(l 



/') -^'^AB 



(98) 



with /' > //3. Thus the noise in the standard form cor- 
responds to white noise, where the noise level (1 — /') is at 
most increased by a factor of 5.5. As in Sec. llV BTI this imme- 
diately follows from the results for the case, where the ideal 
operation is the identity, by simply applying the designed de- 
polarization procedure from the end of the Sec. II V Al with the 
role of particles A' and B' exchanged. Note that the corre- 
sponding twirling procedure remains local w.rt. the physical 
partitioning [AA' ,BB'). 



2. The noisy CNOT-type gate 

We consider now the noisy CNOT-type gate ?7(7r/4), de- 
scribed by the standard form given in Eq. ( I96l l. We will further 
depolarize the corresponding noise process in such a way that 
the fidelity of the ideal operation is decreased (as few as pos- 
sible) and the noise is global white noise, i.e. the simplified 
standard form is given by 

£'p = g[/(7r/4)pt/(^/4)t + (1 - g)l ^ a,a,pa,a, 



qU{n/A)pU{i,/A)^ + 1^1, 



(99) 



where the fidelity of the ideal operation / = g + (1 — 
We find that the amount of noise is increased at most by (ap- 
proximately) an order of magnitude, i.e. (1 — /)/(! — /) ~ 20. 
Clearly, such a further depolarization is only useful if the fi- 
delity of the ideal operation is initially sufficiently large, i.e. 
/ > 0.96, as otherwise the completely depolarizing operation 
would be obtained. 

We will first demonstrate that a depolarization to global 
white noise is possible , and will then discuss the resulting 
decrease of fidelity. The state Es can be written as 



where E is the operator corresponding to the noise part of the 
CPM £s (see Eqs. (|57)) and 

E^/i = |*V4>(*^/4| = J(|*00)(*00| + |*22)(*22| 

+i|*00)(*22| -i|*22)(*00|), (101) 

is the operator corresponding to U{ti/A). 

We will first show that by means of local unitaries, one can 
change the off-diagonal elements of E^j^ in such a way that 
each off-diagonal element in E (or equivalently E) can be 
erased by probabilistically applying either the corresponding 
unitary operation or the identity with appropriate probabil- 
ity. Here, we are no longer restricted to operations that keep 
|5'„/4) invariant, but can use arbitrary local unitaries. We first 

note that by applying Iaa' ® erf (8) erf , one can change the 
sign of the off-diagonal elements |4'oo)(^22| and |^o2)(^'2o|j 
which implies that in the following discussion the sign of the 
off-diagonal elements do not play a role. 

By the depolarization procedure in Sec. IIV B 31 we 
can assume the sub-block Too — spanned by the states 
{|*oo), 1*02), |^'2o), |^'22)}— of the (total) state E to be in 
the form 
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where a = /-too.oo,& = Mo2,02,«u = /ioo,22,iw = /^02.20 (see 
Eqs. (^3). 

We consider the operations Ux = (1 <E) <Jx) and Uz — (1 C?) 
az)- We have that the action of these operations on Bell states 
is given by 



W4|7/-0>,|^l),|V'2),|^3)} 
W{\^0),\^l),\^2),m} 



{\i>i),\iljo),i\i'3),-i\Tp2)}, 



It follows that applying with probability 1/2 the operation 
IaI^x or Ia^^ transforms the state E to the state E^. In 
particular, the subspace Too is transformed to the subspace 
—spanned by {j^-oi), l^-os), |«'2i), |*23>}— and some 
coefficients are aligned. One finds that the resulting elements 
are given by 




where 



^ 01 



C=iA + B)/2, Z=iX + Y)/2, 



(105) 



(106) 



Es = + (1 - f)E, 



(100) 



Note that the element Z is real. At the same time, the subspace 
Foi is transformed to Fqq (where the off diagonal elements are 
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given by i(l — f)w after the transformation, while the diago- 
nal elements are still given by (1 — /)c). Also the subspaces 
Fio and Fn are transformed into each other, where one finds 
that FJq is diagonal with elements (1 — /)e, and also Fjj^ is 
diagonal with elements (1 — f)d. 

Similarly, if one applies the operations W^;, in A instead 
of B, the state E is transformed to a state E"^. In particular, the 
subspace Foo is transformed to the subspace F^q — spanned 
by {|*io), l^-ao), |^'i2), 1*32)}— where the elements of Ffj, 
are the same as of Fqi (see Eq. ( I105» . The transformation of 
the other subspaces follows accordingly, only the role of sys- 
tems A and B is exchanged (e.g. Fqi T^i etc.). Note that 
one can simultaneously change the sign of all off-diagonal el- 
ements of the resulting state E'^ by applying Iaa' ® (8) af . 

If one thus mixes the resulting states E (with probability 
Po), E^ (with probability pi) and E^ (with probability P2), — 
and by appropriately choosing the phases of the corresponding 
off diagonal elements — one can achieve that the final state E" 
has no off diagonal elements in the subspaces Fg\ and F'/q. 
This is guaranteed by choosing 



change the sign of the corresponding off-diagonal elements. 
Choosing 



P=i\X"\ + \Y'- 



(112) 



ensures that the off-diagonal element |^'o2)(^'2o| vanishes, 
while I ^'oo) (4*22 1 becomes 

Y' = ipY" + - p)X" (113) 

and the diagonal elements become 



A' ^pA" + {l-p)B", 
B' =pB" + {l-p)A" . 



(114) 
(115) 



Note that all other elements of E' remain invariant. Finally, 
the remaining off-diagonal element |^'oo)(4'22|, iY', which 
corresponds in part to the ideal operation and in part to the 
noise part, can be formally incorporated into the ideal part of 
the evolution, i.e. the resulting state can be formally rewritten 
as 



Pi = 

P2 = 



z + {\w\ + \x\){i-fy 

\w\il-f) 

z + {\w\ + \x\){i-fy 

1^1(1-/) 
z + (H + N)(i-./)' 



E 



'/=/'l*V4)(*V4l + (l-/')^, 



(116) 



(107) 



The other coefficients of the resulting state can be readily cal- 
culated, taking into account whether a change of sign was nec- 
essary for E^ or E^, where we denote cti = sign(w) + l, 0-2 = 
sign(a::) + 1 with (—1)'^^ = sign(w). One finds that each of 
the subspaces T'^^ , F'/q , r'(^ is diagonal (with all coefficients 
equal) and described by a coefficient, 7oi,7io>7ii respec- 
tively, where 
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po(i - f)c + i-irpic + {~irP2ii - f)e, 

po(l - f)d + (-l)"Vi(l - f)e + {-ir-P2c, 

po(l - f )e + (-l)"Vi(l - f )d + {-Ir-P2{l - 



The subspace Fqq is given by 

PoToo + i-irPiTlo + {-irP2T 



2 

00' 
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where we find that resulting diagonal elements are 

A" = poA + (-l)'^Vi(l - f)c + {-irp2{l - f)d, 

B" = poB + - /)c + {-irP2{l - f)d, (110) 

while the off diagonal elements are given by 

iV" = poiY + - f)w + (-lfHp2{l - f)x, 

iX" = poiX + i-l^Hpiil - f)w + (-l)"^ip2(l - /)a;(lll) 

It remains to show that one can erase the off-diagonal ele- 
ment 1 5*02) (^20 1 > iX. This can be accomplished by randomly 
applying the operation Ia^b or W^ls (see Eq. ( I90» with 
probabilities p and {1 — p). If sign(X) = sign(F), one ap- 
plies in a addition (E) (jf in the second case in order to 



where D is diagonal in the basis {l^'y}} with elements 
dij and /' = 2Y'. The elements djj of D in the blocks 
Foi, Fio, Fii are given by Eq. ( I108> (i.e. doi = rfo3 = c?2i = 
'^23 ~ 7oi etc.), while in the block Fqo we have 



da2 — c?20 — B', 

doo = d22 = A' - Y'. 



(117) 
(118) 



The coefficients dij can even be made equal by further re- 
ducing the fidelity of the ideal operation. In this case, D cor- 
responds to the completely mixed state, and the corresponding 
map is given by global white noise. This is done as follows: 
Using that the probabilistic operation U'^Ib or IaIb, one 
produces a state diagonal in the basis with the same 

diagonal coefficients as Ef, where we consider the situation 
where doo = 0^22 3> dij, i.e. the fidelity of the ideal opera- 
tion is sufficiently large. In this sense, weight from the ideal 
operation can be transferred to the other states. In particu- 
lar, one uses 1^>V^ to increase weight in |^'o2)('I'o2| and 
|'I'2o)(*2o|; ^A^x and ^A^f to increase weight of Fqi; 
U^Ib andUflB to increase weight of Tiq-M^U^ Mfl^f, 
U^Uf and U^U^ to increase weight in Fn. 

One thus ends up with a standard form described by global 
white noise as announced (see Eq. \99\). We will evaluate the 
loss factor for the fidelity for an alternative protocol discussed 
in the next section. This protocol is capable of depolarizing 
also noisy phase gates to a one-parameter standard form. 



3. The noisy phase gate 

In principle, it may be possible to obtain a simplified stan- 
dard form for the gate U {a) with arbitrary a by similar means 
as in the case of the noisy CNOT-type gate J7(7r/4), i.e. by 
manipulating the noisy evolution in such a way that weight is 
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transferred from the ideal evolution to the appropriate noise 
parts. However, for small a one encounters a difficulty which 
may be hard to circumvent. When using parts of the operator 
corresponding to the ideal evolution to eliminate off diagonal 
elements in other parts of the density matrix corresponding to 
the noisy evolution, we have that automatically also the diago- 
nal elements are transferred and hence the noise part is further 
increased. While in the case of U (7r/4), the increase of diag- 
onal elements of the noise part is of the same order of magni- 
tude as the off-diagonal elements, for small a this is no longer 
the case. We have that the off diagonal element of the ideal op- 
erations, Ea, is given by A00.22 = * cos(a) sin(a), while the 
larger diagonal element, Aqo.oo, is given by cos^(a). Imagine 
we have elements in the noise part of order e ^ 1, where both 
diagonal and off diagonal terms of order e appear If one wants 
to eliminate an off-diagonal element in the noise part which 
is of order e, we need (1 — p) cos(a) sin(a) = |e|, i.e. with 
probability (1— p) the off-diagonal element of the ideal evolu- 
tion is transferred to the off-diagonal element of the noise part 
with the sign chosen in such a way that the total off-diagonal 
element in the noise part vanishes. However, by doing such a 
transformation, one of the diagonal elements of the noise part 
is automatically increased by (1 — p) cos^(a) which is of the 
order |e|/ tan(a). Note that for small a, 1/ tan{a) ^ 1, i.e. 
the amount of noise may be increased by orders of magnitude. 
This is clearly not acceptable, as our goal was to obtain a sim- 
plified standard form by sacrificing a relatively small amount 
of the fidelity and not to decrease the fidelity by (several) or- 
ders of magnitude. 



However, under certain conditions one may use an alterna- 
tive method which still allows one to obtain a standard form 
corresponding to global white noise, specified by a single pa- 
rameter In particular, if one can switch the noisy operation 
on and off at will, i.e. one can decide whether one wants to 
apply the noisy operation or does not want to apply it (and 
instead apply something else), then such a depolarization is 
possible. If one considers the case where U (a) is some non- 
local gate, then it is natural to assume such a controllability. 
In this case, one can either apply the noisy evolution with cer- 
tain probability p or apply some other operation with prob- 
ability (1 — p). In particular, one can apply any separable 
operation. This may, however, involve the application of ar- 
bitrary local operations (including measurements), rather than 
the application of local unitaries as we have assumed so far 
Considering the corresponding states, this amounts to mixing 
of the state E with some separable (in the sense A — B) state 
D. The separable operation V associated to D can be im- 
plemented by some random application of local operations, 
V ~ J2i Pi ^iBf piAiBfY . The corresponding Kraus oper- 
ators can be obtained from the spectral decomposition of D as 
indicated in Appendix A. Let us now consider a density matrix 
A of a separable two-qubit state written in the standard basis. 
Then all states of the block diagonal form (see e.g. Eq. (I88> 
with separable block matrices = A for i,j € {0, 1} are 
again separable (recall that Tij denotes subspaces spanned by 
jvPfc;) with k mod 2 = i,l mod 2 = j). In particular, any 



matrix A of the form 




(119) 



with a,(3 G {0, 1, —1, i, —i} as well as any diagonal matrix 
A (with positive coefficients summing up to one) is separable. 
This can be checked by calculating the partial transposition of 
these states, where one finds that the partial transposition is 
positive in all cases which is (for two-qubit states) sufficient 
to ensure separability js^ES- Th^ corresponding separable 
maps can be implemented by a simple sequence of random lo- 
cal unitary operations (in the case of diagonal A), or measure- 
ments (in the case of matrices of the form II 19l . The Kraus 
representation of the state can be obtained as shown in Ap- 
pendix A. 

It is now straightforward to see that mixing E with separa- 
ble operators of the form = A with A given by Eq. 
or an appropriate diagonal matrix, allows one to eliminate all 
(unwanted) off-diagonal elements as well as to align all diag- 
onal coefficients of the noise part. Thus the resulting simpli- 
fied normal form of the noisy operation is given by 



qU{a)pU{a)^ 
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1, 



(120) 



The fidelity of the ideal operation / = q + {1 — q)/16 is 
reduced, where we find e.g. for a = 7r/4 that if / = 1 — e, 
then 



/ > 1 - 17e. 



(121) 



That is, the fidelity of the operation is reduced by about an 
order of magnitude. This can be seen as follows. Consider 
for example the CNOT like gate U{ti/A) with correspond- 
ing standard form of noise E given by Eq. \91\ . We have 
tr(£') = 2a + 26 + 4c + 4d + 4e = (1 - /) and de- 
note y = max (a, 6, c/2, rf/2, e/2), i.e. / < 1 - 2y. One 
can make E diagonal by mixing with matrices F^^ = A. 
Consider for instance the case where all noise is concen- 
trated in 6, V (this in fact turns out to correspond to a (non 
unique) worst case scenario), i.e. y = b. One mixes E 
(with probability p) with a matrix of the form Fqq = A with 
(3 = 0,a = —i (with probability (1 — p)). The resulting 
matrix is diagonal for p — (4w + 1)^^ and has diagonal ele- 
ments b' = bp+ (1 -p)/4 < 2?/(4?; + a' = (1 -p)/4 
and c' ~ d' ^ e' — 0. Note that the worst case corresponds 
to 1) = b. By mixing the resulting matrix (with probability 
q) with a diagonal matrix (with probability 1 — q), one can 
make all diagonal elements equal (which corresponds to white 
noise). We have q > {Ay + l)/(32y + 1). This leads to a total 
final fidelity / = pqf > //(32y + 1) > //(17 - 16/). For 
/ = 1 — e, we thus have 



/>//(17-16/)>l-17e, 



(122) 
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i.e. the fidelity is decreased by about an order of magnitude. 
Note that this formula also holds in the general scenario with 
arbitrary coefficients a, b, c, d, e, u, v, w, x. In the first step, 
the worst case is given when all off diagonal elements are 
maximal, u = a,v = b,w = c,x = d. In the second step 
(making all diagonal elements equal), it is clearly the worst 
case if all weight is concentrated in one element (e.g. b) and 
all others are zero, as one has to fill up the weights of the 
other (14) diagonal elements. Non-zero diagonal elements 
would require less mixing, leading to a larger final fidelity. 
Thus />! — 17(1 — /)isa conservative bound on the final 
fidelity, where in many situations one will end up with a much 
larger final fidelity of the depolarized noisy operation. 



D. Standard form for arbitrary two-qubit unitary operation 

Let us briefly discuss the possibility to generalize the re- 
sults in Sec. IIV Bl and Sec. lIV Cl to the more general case of an 
arbitrary unitary operation as the ideal operation. Consider an 
arbitrary unitary U (or even a class unitaries Ua)- Note that 
similar to Eq. i27\ any unitary two-qubit gate can be repre- 
sented uniquely as ll23ll 

[/^^ = ® Ui e-' S •=! ^''^^^^ (g) (123) 

with J > fJ.1 > IJ.2 > IfJ-sl > and some single-qubit uni- 
taries Ui,Vi {i = 1,2). The main block of the depolariza- 
tion procedures in Sec. IIV Bl was to find a set of {A,A')- 
local unitaries, that leave the corresponding pure state 
(or a class pure states |5'(7q) ) invariant. However, in general, 
— apart from some special cases where e.g. all /i; equal or 
M2 = Ms = — the only local unitary operation that keeps 
the state invariant is given by the identity operation. Hence, 
in generic cases, a depolarization of the operation — under the 
condition that the fidelity of the operation remains invariant — 
is impossible following this approach. Thus no universal stan- 
dard form for arbitrary two-qubit unitary operations can be 
obtained along these lines. 

It seems more appropriate to follow the ideas of Sec. IIV Cl 
and to specifically design such a standard form for a noisy uni- 
tary taking the given form of decoherence into account, and 
allow for an increase of noise. We do not deepen such a dis- 
cussion at this point but rather refer to an alternative approach. 
Instead of regarding the noisy unitary gate as a CPM and al- 
lowing for manipulation before and after the application of 
this operation one might as well consider the dynamical evo- 
lution realizing this gate and allow for a manipulation of this 
evolution by several short intermediate pulses of local unitary 
control operation. Inspired by the results in 138] it is shown 
in the following section that by this procedure it is possible to 
depolarize an arbitrary master equation (of two systems) to a 
standard form described by at most 17 parameters. 



V. STANDARD FORM FOR NOISY EVOLUTIONS 
DESCRIBED BY A MASTER EQUATION 

In this section, we will consider the evolution of two qubits 
described by a master equation of Lindblad form, where the 
ideal evolution is given by an arbitrary two-qubit Hamilto- 
nian. The results can be readily generalized to multi-qubit 
systems whose (noiseless) interaction is described by Hamil- 
tonians that are sums of two-body Hamiltonians. 

We will consider a continuous evolution p{t) = £tp{0) of 
the system due to Markovian quantum dynamics starting at 
t = in the state p(0). Thus the family of quantum operations 
£t forms a Markovian semi-group |l3a|, determined by some 
generator Z, that in the case of two qubits A and B and the 
convention h = 1 satisfies the following differential equation 
{master equation): 

■§^p = Zp := -inp+ -inip + Cp with (124) 
Hp := 

Hip := [H{^^,p] and 

Cp := J2 iki([^i^^P,an + [^i^^,pO). 

k,l 

kT^O V It^O 

We have separated the unitary evolution of the dynamics into 
one part H, which corresponds to the ideal unitary process in 
question, and into a 'Lamb shift' Hi, which is some unitary 
dynamics, that is induced by the coupling between the system 
and its environment and therefore corresponds to noise. The 
more relevant influence of the decoherence process C is how- 
ever incorporated in the 'Liouvillian', which is determined by 
the positive 'GKS' -matrix L = (Lki) 03]. Note that the cor- 
responding sum is over all multi-indices k = (fci,fc2) and 
1 = ili,l2) with = 0, 1, 2, 3 except k (0, 0) or 
1 = 0. Thus the totally mixed state CTq^ = jIab does not 
occur in the sum and L is a positive 15 x 15-matrix. 

Our goal is now to bring some dynamical evolution £t ~ 
e^* into an appropriate standard form £[. Here, £t approxi- 
mates the ideal unitary evolution I given by Ji, where 

Z = -in - iUx + £, (125) 

while the standard form £[ is specified by 

z' = -in-in[ + c'. (126) 

That is, the decoherence process corresponding to the stan- 
dard form is described by Hi and L', rather than H\ and L in 
the original evolution. As in the case of CPMs, our goal is to 
obtain a simplified standards form in the sense that the num- 
ber of relevant parameters describing the decoherence process 
are decreased, while the desired Hamiltonian evolution is not 
altered. To achieve this we will make use of the following 
facts iH: 

(i) Let U be some unitary matrix and Up := UpW be 
the corresponding operation. By unitary conjugation 
the Markovian evolution £t = e^* can be transformed 
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into the Markovian evolution £[ = U o £t o U'^ = * 
described by 

H' = UHU^ , (127) 

Hi = UHU^ , (128) 

L' ^ OLO^ , (129) 

where O is the orthogonal matrix corresponding to U, 
that describes the basis change tTk i-^ Ua-^JJ'^ for the 
linear basis of hermitian traceless operators dk- 

(ii) A Markovian evolution * according to a lin- 
ear combination Z' = Yl,^^iPi^i (EjPi 1) 
can be simulated by repeatedly applying a sequence 
e-2iAt . . . g^nAt gojj^g small time intervals At = jj 
(M is the number of repetitions), i.e. 

(ntieP-^-^)*' eSf=iP.2.t. (130) 

Note that the approximated GKS matrix is L' = 
"^■PiLii and the error in approximation is of order 
0(Ai2) 0|. 

An alternative method consists in the random applica- 
tion of the time evolutions e^'^* with probability pi in 
each of the time intervals At = t/M. That is, the evo- 
lution in the time interval At is given by X^i-P*^^'^*' 
which accurately approximates the desired operation in 
first order At. A sequential application of these ran- 
dom operations M times reproduces the desired opera- 
tion i*=i P'^'* in the Hmit M oo, i.e. 0| 

{EtiP^e^'^y ^t^va.t . (131) 

Following the structure of the previous section we first con- 
sider the case of decoherence itself, i.e. we set H — 0. 
We propose a depolarization protocol that, after integration 
of the corresponding master equation, yields the same stan- 
dard forms as obtained in Sec.|lll]and Sec. lIV Al Second, we 
consider the case where H corresponds to some Ising-type in- 
teraction, e.g. H = (jy ® (jy . We make use of the results 
obtained in Sec. lIVB"2l for the corresponding unitary e~'^*. 
We show that a depolarization procedure exists for which — 
in the limit of infinite intermediate local control operations — 
the system evolves according to some standard form which 
has the standard form Eq. ( I88> . when regarded as a CPM £[. 
The depolarization procedures we describe in the following 
(Sec. IV Al and Sec. IV Bt are only relevant in cases where one 
is interested in the standard form for the complete dynamics 
and not only after some given time (the latter corresponding 
to the case of CPMs discussed in the previous Sections). Oth- 
erwise, one may use the conceptually simpler depolarization 
procedure for CPMs. Nevertheless subsections IV Al and IV Bl 
provide the necessary tools for the procedure proposed in the 
subsequent Sec. IV C| where we show how to achieve a stan- 
dard form for some arbitrary unitary dynamics. Although 
this procedure overcomes the restrictions to the area of ap- 
plications in Sec. IIVBI this depolarization protocol generally 



increases the noise level of the decoherence process. In the 
following we assume that local unitary control operations can 
be performed on time scales negligible compared to the in- 
teraction time for the dynamics. We will thus refer to these 
operations as being instantaneous. 

A. Standard form for decoherence processes 

We first consider maps describing pure decoherence pro- 
cesses {H = 0) for a single qubit. For the depolarization we 
consider the same twirling procedures as in Sec.|n] but now 
we intend to bring the Markovian generator Z of the initial 
dynamics into the standard form 

Z' = ^ UkZUl , (132) 

k 

where Uk denote the unitaries which were applied in Sec. 
with equal probability Uk to achieve the standard form of a 
Pauli channel, i.e. Uk = cTk is one of the Pauli matrices and 
Uk = \, or the simpler standard form of a depolarizing chan- 
nel, i.e. Uk is of the form Qiai (Qi = e'^'^' , Z = 1, 2, 3) and 
Uk = In Sec.|III|the number Uk represented the probabil- 
ity to apply the different twirling unitaries. According to (ii), 
a similar procedure still works, where intermediate random 
applications of the corresponding unitaries lead to a standard 
form of the Markovian generator Z. Alternatively, one can 
consider a further splitting of the time interval Deltat (see 
Eq. I130> . where all possible unitary operations corresponding 
to the depolarization process are applied sequentially. We will 
consider the second approach in the following. 

More precisely, we will consider the following depolariza- 
tion protocol: Let the actual dynamics of the system (i.e. the 
decoherence process) evolve for some time t and choose a 
split of the total time t into M sufficiently small time intervals 
At. During each of these small time intervals the system dy- 
namics is accompanied by the sequence of instantaneous local 
control operations Uk+iUk (Uq = 1) applied in arbitrary or- 
der but with equal distance Uk At. From fact (i) it follows 
that instead of the original dynamics e*^ ^* during each of 
the time intervals Uk At the system evolves according to the 
Markovian generator UkZU\. In the time interval At the evo- 
lution is thus given by 

-Q g«fcAtw,2wt ^ (J33) 

fe 

If these intervals are chosen sufficiently small (i.e. A4 —* 
oo) fact (ii) implies that the overall system dynamics can 
effectively be approximated by the Markovian generator in 
Eq. (fT32li . 

Let us consider the effect of this depolarization on L and 
H\ more closely. According to (i) the GKS matrix is brought 
into the standard form Eq. OOt or Eq. i35i . except that the first 
row and column of Eq. ( I29> is disregarded in both equations. 
Thus we obtain 

3 

£'p^^2Y,Lk{p~<ykP<Jk) (134) 

k=l 
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with Lfe = Lkk in the case of twiriing with PauH operators 
Uk= Ok and 



C'p 



-2L 



o-kpo-k 



-AL (2p- 1) (135) 



fe=i 



with L = i (ill + L22 + -^^33) in the case of complete de- 
polarization. Similarly for the Hamiltonian H\ the twirling by 
means of the Pauli matrices Uk — <Jk yields (see (i)): 



1 ^ 1 
: - <^Mak = - tr{Hi) 

k=0 



(136) 



Since this twirling is also performed before each of the Clif- 
ford unitaries Qi applied, we obtain the same result in the 
case of complete depolarization. Thus in both cases the Lamb 
shift Hamiltonian in the standard form gives only rise to some 
overall phase factor e~5"(^i) which can be neglected. 

We briefly examine the dynamics E'^ due to the Standard 
forms, i.e. the solutions of the master equation 



p = C'p , 



(137) 



where the Liouvillian C is given by Eq. or Eq. ( I135> . 

It is straightforward to see that in case of Eq. ( I134t the Pauli 
matrices diagonalize the Liouvillian C |42], i.e. 

C'a^ = C'ai = -4(^2 + ^3)^1 

£'0-2 = -4(Li + L3)(T2 C'az = -4(ii + ^2)^3 • 

For an arbitrary initial state p(0) = ^ [1 + n • ct] (|n| = 1, 
a = (i7i, (72, (Ta)"^) we thus obtain 



withp(i) = e"^^*. 

The fact that the above depolarization procedure for contin- 
uous dynamics rediscover the standard forms already obtained 
in Sec.|lIl]for CPM can also be understood by the simple fact, 
that the evolution * is invariant under unitary conjugation, 
i.e. Ue^ ^U^ = *, iff the generator Z' is invariant under 
this conjugation. By construction of the protocol above, this 
definitely holds for the standard form Z' w.rt. any of the cor- 
responding twirling unitaries Uk- With Uk being the twirling 
constant we thus obtain 



E' 



„z't 



k ■ 



(141) 



Since this Eq. precisely reflects the corresponding depolariza- 
tion procedure for CPM of Sec.|in| the resulting time depen- 
dent CPM £^ has to be of standard form for CPMs. 

The depolarization protocols for decoherence processes 
(described in this subsection) can be readily generalized to 
the multi-party setting. The Lamb shift HI in the standard 
form can be neglected again whereas the Liouvillian is of a 
standard form that corresponds to the multi-party Pauli chan- 
nel Eq. dsn or the multi-party depolarizing channel Eq. ( I53l l. 
Note, however, that each additional party causes a finer split- 
ting of the time interval At, yielding an exponential increase 
of the number of intermediate control operations. More pre- 
cisely for N parties the system dynamics in each time interval 
A has to be interrupted by 4^ [12^] control operations in or- 
der to achieve a dynamics corresponding to a Pauli channel 
[depolarizing channel]. In this case, the alternative method of 
applying random local unitary operations rather than the com- 
plete sequence of unitaries is certainly privileged. 



p(t)=e^'V(0) = ^[l + n(t).a] 



1 



(138) 



where n(t) = (nie-''(-^^+^-^)*, 7126-4(^1+^3)*^ „gg-4(Li+L2)t 
The action of E'^ in terms of Pauli matrices as in Eq. ( I29t is 
actually that of a Pauli channel 



B. Standard forms for noisy Ising-type interactions 



^ Let us now move to the case where the ideal operation is 
I given by some Ising-type interaction 



H 



AB „^AB 



(142) 



where 



Etp^Y Ek{t)cTkpCTk 



(139) 



fe=0 



Mt) = 4 (l + e 
Mt) = 7 fi + e 



4(L2+i3)t I „-4(ii+-t3)t I ^-4{il+i2)t 



+ e 



+ e 



-4(L2+i3)t „-4(ii+i3)t ^-4{-tl+i2)i 



E2{t) = i ^1 _ e-*'^2+i3)t ^ g-4(ii+i3)t _ g-4{Li+i2)t 



T 1-e 



-4(L2+i3)t _ g-4(il+i3)t _|_ g-4{ii+i2)t 



Similarly a specialization to complete depolarization reveals 
the depolarizing channel 

£t{p)=p{t)p+{l-p{t))\x{p)]^l (140) 



For depolarization we can essentially consider the same pro- 
tocol as in the previous section, except that we now take the 
twirling unitaries Uk, that were used in Sec. lIV B"^ More pre- 
cisely Uk ~ and the Uk are given by the 32 unitaries from 
the product set U\-U2- U^, where 



y , e 



/I 1 ^^1 1 ABx 

V-Ai-B, CTy i-B, ^A<^y , j , 



- 

V e 



} 

(143) 



e.g. Uk = e~''T°'yax <E) (TxCry- Recall that a twirling with 
these unitaries brings the phase gate U (a) = e"*^" with an 
arbitrary angle a into a standard form Eq. ( I88> described by 
only 17 independent parameters. 

Apart from the different choice of twirling unitaries the de- 
polarization reads exactly as in Sec. IV Al The overall inter- 
action time t is divided into sufficiently small time intervals 
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At, in which the system dynamics is interrupted for short lo- 
cal unitary control operations Uk- The resulting dynamics ap- 
proximates an evolution £^ — * with new Markovian gen- 
erator 

32 

Z' ^Y.''kUkZUl^-tH-tH{ + C' . (144) 

Note that due to the choice of the twirling unitaries the ideal 
Ising-type interaction Hamiltonian H' = H is in fact not 
changed by this protocol. The GKS matrix L' is of standard 
form Eq. ( I88t . except that in Eq. J88> the first row and col- 
umn are disregarded. Since no normalization constraints are 
involved, the smaller 15 x 15-matrix is still specified by 17 
real parameters. For the new Hamiltonian of the lamb shift 
one can compute that the twirling yields H( = Hl^QaQcr^ + 
H\2<^2'^2 ■ Omitting the term T/qq Iab, which would con- 
tribute only as global phase 6^*^°°* to the system dynamics, 
the Lamb shift in the standard form thus can be specified by 
a single real parameter H\2- Moreover we will include this 
term into the ideal interaction term rewriting H' ~ g'ayCTy 
with g' = g + i?22 ^iid will disregard any Lamb shift in the 
following. The change in coupling strength demands a rein- 
terpretation of the interaction time. In order to simulate the 
actual dynamics e^* running for some time t by a dynamics 
of the standard form * the simulation actually has to run 
for the time tg = ct, where c = ^ represents the time cost of 
the simulation. 

A similar argument as in the previous subsection shows that 
the obtained standard form * seen as a CPM is actually in 
the standard form considered in Sec. IIV B 21 



C. Standard forms for arbitrary noisy evolutions by means of 
sacrificing 

Let us now consider standard forms for arbitrary ideal uni- 
tary evolutions H. We make use of the fundamental fact 13^1 
that by a stroboscopic application of a sequence of local uni- 
taries any (entangling) two-qubit Hamiltonian H can simu- 
late the Hamiltonian Hy ~ ay ® <jy of the phase gate oper- 
ation U{a) = e"'^"" in Sec. IIVB2l to arbitrary good ap- 
proximation (and vice versa). Before going into detail we 
shortly sketch the procedure of deriving a standard form for 
any noisy unitary evolution. Along each single stroboscopic 
step At = jj {M number of steps) of the simulation protocols 
we propose to 

(1) first apply the sequence of unitary operations in order to 
obtain the evolution according to —iHy — iHiy + Cy, 

(2) depolarize the CPM Tt described by —iHy — iH\y + Cy 
according to the procedure derived in Sec. IV Bl vielding 
a standard form !F[ given by —iHy — iHly + Cy and 
finally 

(3) transform it back to the original Hamiltonian H 
accompanied by some decoherence process of desired 



standard form 7Y( and C . 



approximates simulate 



standard form standard form 

n < z' < z' 

approximates simulate 

We remark that steps (1) and (3) require in general a 'time 
cost', i.e. the simulation of the action of a desired Hamilto- 
nian for some time ct requires a time t. This time cost trans- 
lates into a smaller pre-factor for the interaction Hamiltonian 
in the corresponding generator Z', ultimately leading to an 
increased noise. That is, the ratio of the strength of desired in- 
teraction (described by TC) to strength of noise (described by 
C) decreases, leading to a reduction of the fidelity. We note 
that for two-qubit systems the time cost is at most 3. 

In order to simulate the Hamiltonian Hy by the (entangling) 
two-qubit Hamiltonian H and fast local unitary transforma- 
tions (see 13^ ') one considers the decomposition of Hy in 
terms of H (term isolation): 

1=1 

where Vi are the local unitaries with probabilities Vi > 
Vi — 1), Qi and Q2 are some local Hamiltonian on qubit 
A and B respectively and > is some factor to adjust 
the coupling 'strength' of the Hamiltonians H and Hy. If 
the unitary evolution e~*^f * is supposed to be simulated for 
the time t, the simulation has to be carried out for the time 
tg = Cyt. Since the local unitary control operation can be 
performed on negligible time scales, the factor Cy thus deter- 
mines the time cost for the following simulation: One chooses 
a split of the time tg into M time intervals At, such that the 

sequence e-™'^''«^'^* ^-'QiAt g-^Q^At j ^ ^^f^. 

cient approximation for e~*^!'*= as discussed in (ii). Note that 
in each time step At the original dynamics according to H is 
simply interrupted after the time Vi At in order to apply the lo- 
cal unitaries V^!|_i (Vq = 1). This corresponds to the system 
evolving according to a sequence of Hamiltonians ViHV^ for 
the time intervals ViAt. At the end of each simulation step At 
the local unitary e~*2iAt ^ g-jQ2At j^^^ j.^ applied in or- 
der to cope with the single qubit dynamics in the Hamiltonian 
simulation. 

Similarly, one can consider a Hamiltonian simulation for 
step (3) according to the decomposition 

with local unitaries Wj, single qubit Hamiltonians Qi, Q2 
and time cost Cm for the 'backward' simulation. For step (2) 
we use the twirling protocol derived in Sec. IV Bi providing a 
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standard form Eq. (I88> described by 17 independent parame- 
ters. With Cy we denote the corresponding time cost of this 
depolarization procedure. 

With these notations at hand we can now specify the proto- 
col to achieve the standard form E'^ — * for an arbitrary 
noisy two-qubit evolution £t = e^*. Let £t = e^* be a 
Markovian evolution with the generator Z = —iH — iTi.\ + C, 
where H p = —i[H, p] corresponds to the ideal evolution with 
Hamiltonian H,T-C\p = —i[H\,p] represents the Lamb shift 
with Hamiltonian H\ and 

Cp = ^ iki ([c^kP, cTi] + [(Tk, pcri]) (147) 

corresponds to the Liouvillian with GKS-matrix L. For nota- 
tional simplicity we will in the following restrict to the case 
where in both steps (1) and (3) no single qubit dynamics has 
to be corrected, i.e. the terms Qi, Q2, Qi and Q2 in the de- 
compositions Eq. ( II45I 1 and Eq. ( I146> vanish. If the system 
evolves for some time t, the following protocol requires the 
time ts := CyCyCwt and thus has time cost c^CyC^. The 'sim- 
ulation' tg again has to be divided into sufficiently small time 
steps At = |j. In these time intervals we consider the fol- 
lowing sequence of i? = 32 Ru, operations: 

R 

Y[ WjUkV^e^'"'''''''^^'vlulw] . (148) 

i.j,k—l 



As discussed in Sec. IV BI the effect of step (2) on Lj, is to 
bring the matrix into the standard form Eq. (I88> . The final 
standard form L' of the GKS matrix is obtained from by 
mixing according to {wj, Wj) and is thus specified by 17 in- 
dependent parameters only, although L' in general is not of 
the form Eq. ( I88> . As seen in Sec. IV Bl we can neglect the 
lamb shift by introducing some time cost Cy. 

To summarize, we have shown how to achieve a standard 
form for arbitrary two-qubit interactions, where the noise pro- 
cess (described by the GKS matrix) is specified by 17 param- 
eters. The above protocol can be affected by different sources 
of errors. For this, one can again compare the noise level 
of the standard form dynamics Z' with the noise level of the 
original dynamics Z in terms of the distance d{£t,It) and 
d{£f,Tt) to the ideal unitary evolution Tt = e*'^* for differ- 
ent times t, where d{£,T) is a suitable distance measure (see 
Sec. IIIO . Although we have yet not performed a detailed 
error analysis in this sense, a non-unit time cost (at most a 
factor of CyCw < 3 from simulating corresponding Hamilto- 
nians - steps (1) and (3), plus the time cost Cy from 'Lamb' 



This sequence of operations corresponds to a splitting of the 
time interval At into smaller intervals of length WjUkVi At, 
in which at the beginning the (fast) local unitary WjUlV^ 
is performed, the system then evolves according to the given 
dynamics Z and finally the inverse unitary 'pulse' WjUkVi 
is applied at the end of the interval WjUkVi At. In the limit 
At we obtain the Markovian dynamics £[ = * with 



Z' = WjUkV,W,UkV.,ZVlulw] . (149) 

ij^k—l 



It is straightforward to show that the ideal operation Ti, in the 
generator Z' remains the same, since the twirling over Uk 
leaves the Hamiltonian Hy invariant. Moreover Z' again has 
a decomposition of the form Z' = —iH — iT-C[ + C described 
by the new lamb shift Hamiltonian HI and the GKS matrix 
L', that are obtained as follows: 



Z' 

= Y.,^jW,H'^,WU] (150) 



shift), in general, corresponds to an increase of the noise level 
for the evolution. 

As it holds for the depolarization of CPMs in previous chap- 
ters and as opposed to the assumptions made in this paper, in 
practise, the depolarization protocol has to face imperfections 
in the local control operations, whose extent depends on the 
physical realization. Additionally, for the depolarization of 
master equations by means of stroboscopic control operations 
one also encounters errors of order O(At^) due to the finite 
approximation (see fact (ii)) j40ll . Note that, in practise, there 
will be a trade-of between errors in approximation and errors 
due to imperfect local control operations. 

D. Simplified standard forms for arbitrary noisy evolutions 

A further reduction of the number of relevant noise param- 
eters and thus a simpler standard form may be achieved fol- 
lowing the ideas developed in Sec. lIV Cl for CPMs. There, by 
increasing the noise level and hence reducing the fidelity of 



Z 



Vi.Vi 



(2) 



5=" 



I 



(3) 
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the operation, we have shown that one can in fact achieve that 
the noise part of the evolution is described by only a single pa- 
rameter (white noise). The procedure outlined in Sec. live 31 
is based on probabilistically mixing the (already depolarized) 
noisy CPM £ with a certain separable map T), i.e. a map 
which can be obtained without interactions between particles. 
That is, one chooses randomly whether one wishes to apply 
the map £ corresponding to the noisy operation, or the sepa- 
rable map T). For a proper choice of T) the resulting map £s 
is of the form Eq. ( I120> . 

In the case of master equations, one may adopt this proce- 
dure in such a way that for each time interval At, one applies 
the (already depolarized) noisy evolution described by the Z', 
together with an appropriate separable evolution (that may e.g. 
be generated with the help of available local unitary control 
operations and additional measurements) with corresponding 
generator Zd- Both evolutions now have to be apphed ei- 
ther sequentially or chosen randomly. This implies that either 
one has the ability to switch off the evolution Z', or one can 
produce a separable evolution of arbitrary strength Zj^. Note 
that in this case, fast local unitary operations are in general 
not sufficient, but arbitrary local control operations (including 
measurements) are required to generate the desired separable 
operations. As in the case of CPMs this depolarization proce- 
dure requires moreover the knowledge of the exact form for 
Z' in order to choose an appropriate, separable Zo- The total 
evolution is finally described by a Liouvillian with GKS ma- 
trix proportional to the identity, which corresponds to global 
white noise at the level of the respective CPM. 



We have generalized our results to evolutions described by 
a master equation of Lindblad form. Standard forms for deco- 
herence processes and interaction Hamiltonians proportional 
to the Ising Hamiltonian can be derived using similar methods 
as for CPMs, leading to standard forms with same number of 
parameters. We have also obtained a standard form described 
by 17 parameters for arbitrary two-qubit interaction Hamil- 
tonians, which, in general, goes along with an increase of the 
noise level. As the basic tool we have used the possibility 
to simulate the Ising Hamiltonian by an arbitrary Hamiltonian 
(and vice versa), together with depolarization of the Ising type 
interaction. Again, a further simplification to a single parame- 
ter leading to a GKS matrix proportional to identity is possible 
under certain circumstances. 

We are confident that such simplified standard forms for 
noise processes will provide a useful tool to investigate var- 
ious problems in quantum information processing involving 
noisy apparatus and interactions with environment. Straight- 
forward applications include the possibility to calculate lower 
bounds on the channel capacity of arbitrary noise channels (by 
investigating the corresponding depolarized channels), and a 
simplified process tomography where only a reduced number 
of parameters of the noise process needs to be determined. 
Further conceivable applications include the determination of 
lower bounds to the lifetime of entangled states, and strict er- 
ror thresholds for quantum computation that are valid for ar- 
bitrary noise processes and are not restricted to certain noise 
models. 



VI. SUMMARY 

In this article, we have introduced the concept of depolar- 
ization of noisy evolutions. We have shown how to reduce 
the relevant number of parameters describing an arbitrary, un- 
known noise process described by a CPM in such a way that 
the ideal (unitary) part of the evolution is not altered. For de- 
coherence processes we have explicitly calculated the corre- 
sponding standard forms for multipartite systems of arbitrary 
number N of parties and arbitrary dimension d. We find a re- 
duction of an arbitrary noise process described by 0{d'^^) to 
local and global white noise processes described by only 2^ 
parameters. For specific two-qubit unitary operations (e.g. 
Phase gate with arbitrary phase), we obtain a standard form 
described by at most 17 parameters. For other gates, the stan- 
dard forms can be further simplified. In particular we find 
standard forms described by 8 parameters for the CNOT-gate 
and 3 parameters for the SWAP gate. The depolarization pro- 
cedures used to obtain these standard forms are universal in 
the sense that the exact form of the noise process need not 
be known. With knowledge of the exact form of the noise 
process, and by allowing for a (small) reduction of the fi- 
delity, one can further simplify the standard forms. In fact, 
we have derived a depolarization protocol that yields a reduc- 
tion to global white noise, which is described by only a single 
noise parameter and where, in the worst case, the noise level 
is increased by about an order of magnitude. 
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Appendix A: Spectral decomposition versus Kraus 
representation 

For sake of completeness we will shortly discuss the rela- 
tion of the Jamiolkowski state E with two common represen- 
tations of the corresponding CPM £ in terms of its Kraus rep- 
resentation and its purification. In Appendix A we consider a 
CPM £ described by its Kraus representation. More precisely 
any CPM £ acting on H"^ allows for a decomposition of the 
form 

r 

£{M) - K.MKI (151) 

i=l 

with r < dA dA' (r = Choi rank) Kraus operators Ki € 
A^(H"^, H"* ), where Ki can be chosen to be orthogonal, i.e. 
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5ij . £ corresponds to SLOCC operation iff 



K^K, < 1 



(152) 



where equality holds iff £ is trace preserving. Eq. (I151> is 
an immediate consequence of the corresponding fact, that any 
positive operator E' > (see No. 3 in Sec. Ill B> allows a spec- 
tral decomposition 



Appendix B: Purification for quantum states and operations 

Frequently a CPM £ is also regarded as a description of the 
(non-unitary) evolution of the system alone (i.e. by tracing the 
environmental degrees of freedom), where the system A to- 
gether with its environment C is believed to evolve according 
to some unitary operation Ug . Any such system-environment 
model Ug then is said to be a purification of the CPM £, if it 
yields £ as the evolution of the system alone after tracing out 
the environment: 



E = 



(153) 



£{p)=tvc [U, 



AC 



PC (U, 



Ac^r 



(158) 



where IvA are some unnormalized vectors in H 



A' , 



rank(i?) < x dA' , which can be chosen to be orthogonal. 
Using the decomposition 



E 



ct/3| \A' 



this correspondence is simply given by 

K , 



a/3 I ^ A' A, 



or K^ ' = v«A v.^ . 

Tracing out either system A' or A yields: 



(154) 



(155) 



Here we assume that the system A and its environment A' 
are initially decoupled p (g) pc, where pc = is the 

maximally mixed state of the environment. Like the op- 
erator sum decomposition, a purification of a quantum op- 
eration £ in terms of unitary evolution Us can simply be 
derived from the corresponding purification of the quantum 
state E in terms of a pure state on the joint system A' A 
and C. For this let \ip)^ ^'^ be the pure state, such that 



E 



A' A 



, then the corresponding purifi- 



cation JJf is given similarly to Eq. ( I155> by 



C//^ = ^/dA X dc 



ijk 



V A' A 



{j\^{k\, (159) 



K, 



^A'E^'^ = y: [ki 

i=l 



(156) 
(157) 



using the decomposition I^A) = Z^i.j.fc 

Note that in both cases purifications can be chosen such that 

dc = dimcili'-^) < dA x dA'- 



Note that No. 4 in Sec. IllBI now follows directly from 
Eq. ( I152t using Eq. ( I156> . The unitary freedom in the choice 
of decompositions for E translates to the corresponding CPM 
as follows. Two decomposition of £ [E] with Kraus opera- 
tors Ki and Lj [vectors \vi) and \wj)] correspond to the same 
CPM [positive operator] iff there exists a unitary matrix Uij 
such that Ki = J^ij UijLj [\vi) = Y,ij Uij \wj)]. From these 
decompositions one can easily obtain a few more results about 
the relation between positive operators and CPM under the 
Jamiolkowski isomorphism: 

10. £ is facto rizable, i.e. £{M) = KMK\ iff E is pure, 
i.e. E = 

11. £ is an isometry, i.e. £{M) = VMV^ with V^V = 
\a, iff dA < dA' and E is maximally entangled, i.e. 
pure and Ma'E^ = 1a. (For a unitary we have dA = 

dA'-) 

12. £ is a projection, i.e. £{AI) — dA tr [t|M] ti, iff E is 
a product state, i.e. E ^ ti ® T2- 

13. £ can be decomposed into a sum of projections, iff E is 
separable. In this case £ is entanglement breaking j43ll . 



Appendix C: Details for the extension to d-level systems 

In this appendix we present some details about the general- 
ization of the standard forms for decoherence to d-level sys- 
tems as it was introduced in Sec. |lll] We proof that 

(i) the twirling over the Pauli group as in Eq. (I41> depo- 
larizes any CPM to the standard form of a generalized 
Pauli channel Eq. (I42> . 

(ii) the standard form of a generalized depolarizing channel 
Eq. ( I50> can be achieved by twirUng over a finite set of 
generalized Clifford operations. 

But let us briefly consider the measurements in the gen- 
eralized Bell basis Eq. ( I36> of the isomorphism protocol. It 
is straight forward to compute, that the other Bell measure- 
ment results \il}ki)'^^ yield £ (U^ipU'^i) instead of £{p) for 
the outcome of the protocol on system A' but with the same 
probability Jj. 

Before coming to the proofs, we summarize some useful 
facts (see also ll27ll ) about the generalized Pauli operators Uki 
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(see Eq. i37\ ). that are straightforward to proof: 

Uk'i'Uki = e'~*"' ■'C/(fc+fc')(;+;') 



■'(fe)(-0 
e'^'''''lV'(/c+fc')(/+;')> 



-,;2zL(fc+fc')-r 



k+k')(l-l')l 



(160) 

(161) 
(162) 

(163) 

(164) 

(165) 



With these relations at hand it is easy to verify (i). 
First note that the set of generalized local Pauli operators 

S = {gki\ 9ki = ® U^lkm ; ^ e } (166) 

is a commutative subgroup (of the generalized local Pauli 
group) that stabilizes |0) and is generated by the two elements 
giQ and goi- Moreover a simple calculation shows that 

Sfe'rlV-fei) =e'^(''''"'''')|^fei> • (167) 

For a general state E I]q/3,q'/3' Eafj.a'p' iV'a/j) (V'a'/J; I we 
therefore find, that it can be diagonalized by a probabilistic 

application of the local unitaries g^i with uniform probability 
1 . 



^ d-l 



l2 E Eal3,a'0' \'4^ap){ll^a'l3'\ 



X 



ap.a'fj' 



X I > e I ^ 

\k=o I \i=o 

= E ^a/3,a/3 I V'a/3> ('/'a/3 I 

Q/3 



(168) 



This mixing operation V — Vi o 'D2 can also be decomposed 
into a mixing of shift operation 'Di{E) = i Ylfl^ goiEg^i 
and a mixing of phase multiplication operation 'D2{E) = 

Let us now consider statement (ii). 
For this we assume that £ is already brought into the form of 
a (generalized) Pauli channel of Eq. ( I42> by a random appli- 
cation of one of the Pauli operators Uki- Now the generalized 
Clifford operation is a unitary operation Q that maps the gen- 

eraUzed PauH group V = ^e^^^Uki\k,l,S €Nd^ to itself 

under conjugation, i.e. QVQ^ ~ V. It is totally specified 
by its action on the two generators Uw and Uai of the Pauli 
group, since 



QUkiQ^ = {QUoiQ^y (QUioQ^)' 



(ka+lb)(kc+ld) ■ 



(169) 
(170) 



Here we have chosen QUiqQ^ = Uac and QC/oiQ^ = Um- 
In addition we disregard appropriate phase factors e*^^ with 



5 S N2(j, since they will be irrelevant for our purposes. Up 
to these phase factors the Clifford unitary Q permutes a Pauli 
operator Uki to a new element Uk'v, which in modular arith- 
metic (modulo d) is related to Uki by the linear transformation 



a b 
c d 



(171) 
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This linear transformation needs to be sym plec tic 14411 in order 
to truly correspond to a Clifford operation 14511 . Symplecticity 
of Cq in our single-party case simply reduces to the condition 
detCg ~ 1 (modulo d). Applying one of these CUfford uni- 
taries Q"^ (g) Q*^ to the state E^"^' gives the state 

= E Eki\^k'i')^^'{^k'i'\ , (172) 

where as claimed above any phase factor would cancel out. 
Note that the component Eqq will remain 'untouched', since 
any symplectic matrix Cq (even over F^; with d non-prime) is 
invertible. In the following the CUfford unitary Q will be cho- 
sen uniformly at random from a set of all Clifford unitaries, 
where each Q corresponds only to a single Cq (i.e.: fix a 
choice of phase factor for each Cg). 

By elementary results from group theory it follows, that by 
application of the different C the set of all vectors (fc, l)'^ with 
/c 7^ or / 7^ will be mapped onto itself in such a way, 
that all vectors will occur equally often. For this let G de- 
note the group of symplectic matrices over F^, that act on the 
set X — F"^. Furthermore for x G X let Gx := {gx\g € 
G} denote the orbit of x under the group action G and let 
Gx ■= {g & G \ gx = x} denote the stabilizer of x. From the 
stabilizer orbit and Lagrange theorem it follows that for any 
finite group G acting on a set X we have \G\ = \Gx \ \Gx\. 
This result can be used to show that each non-trivial element 
y e X\{Q}iQ := (00)^ G F^) is obtained \G\ times, if the 
complete group G is applied to all elements in X \ {0}. Since 
G consists of invertible matrices, it maps the set X \ {0} onto 
itself. Moreover any y G X\{0}is only obtained from el- 
ements of its orbit Gy. For a fixed element x £ Gy the set 
Gxy :^ {g G G \ gx = y} can be rewritten in terms of only 
one of its elements g' (i.e. g'x ~ y) and the stabilizer Gx as 
Gxy = g'Gx- Thus y is obtained form each of the \Gy \ ele- 
ments (of its orbit) by \Gx \ different matrices. Since for two 
elements in the same orbit we have Gx = Gy, any element 
y is obtained \Gx \ \Gx\ — \G\ times. Note that in the case of 
prime dimension d, the set X is (not only a module but also) a 
vector space over the field F^; and one can easily show that for 
all a; 7^ the orbits are the same Gx = ^ \ {0}. This is due 
to the fact, that for each non-trivial vector x one can find a 
symplectic (i.e. F^-invertible) matrix g with the first column 
being x (and the second the orthonormal vector x^). 

A random application of the corresponding Clifford opera- 
tions therefore provides a mixing of all components Eki with 
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^ or / 7^ 0. Thus starting with a CPM £ in the form of 
a Pauli channel ( Eq. J42l i) we can achieve the standard form 
in Eq. ( 15 Ot by uniformly choosing a unitary Q form the set 
of Clifford operations and applying before and Q after the 
application of the CPM £. In fact the actual set, which the 
Clifford operations have to be chosen from in order to achieve 



a complete mixing of all the components E^i with k ^ Q or 
I 7^ 0, might even be decreased, as it is illustrated for the qubit 
case in Eq. ( I33> . Note that the Clifford operation Qi, Q2 and 
in Eq. ( I33> correspond to the three symplectic matrices 
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